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We describe a method to derive, from first principles, the long-distance asymp- 
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1 Introduction 



The aim of this article is to derive the long-distance asymptotic behavior of correlation functions 
of integrable models "from first principles". To reach this goal, we will use the framework of the 
Bethe ansatz [Jl [75l [80l [Ml [87] in its algebraic version |[28l [79] . For simplicity, we shall explain 
the main features of our approach using the example of the XXZ spin-^ Heisenberg model [36117]. 
However, it will become apparent that the method we propose here is quite general and applies, 
for example, even to continuous integrable models like the interacting one-dimensional Bose gas 

The starting point of our approach is an exact expression, that we called the master equation 
representation, for the generating function of the two-point correlation functions of these models 
that we obtained in our previous works [IH [50l [5TJ [l5] . In the case of the XXZ spin-^ Heisenberg 
model (at non zero magnetic field and in the disordered regime), this means that we will begin 
our investigation from the exact expressions for the correlation functions on a finite lattice, and 
then take successively the thermodynamic limit and the large distance limit. The long-distance 
asymptotic behavior that we obtain confirms, at the leading order, the predictions given from 
bosonization [67l [33l [Ml [35] and conformal field theories [TTl fTl fTHl [8] . 

1.1 Historical context: a brief survey 

Bethe ansatz and especially its algebraic versions [7llZ5l|8ni[E6l[87l[28l[Z9l|6l[29l[63l[5^ provide a 
powerful framework for the construction and the resolution of a wide class of quantum integrable 
models in low dimension. The central object of this approach is the i?-matrix satisfying the (cu- 
bic) Yang-Baxter equation and providing the structure constants of the associated (quadratic) 
Yang-Baxter algebra for the operator entries of the quantum monodromy matrix T{X). The 
representation theory for this algebra leads to the construction of integrable models operator al- 
gebras including in a natural way the Hamiltonian, its associated commuting conserved charges 
generated by the transfer matrix T(A) = trT(A) together with the creation-annihilation opera- 
tors determining their common spectrum. 

The archetype of such models is provided by the XXZ Heisenberg spin-^ chain in a magnetic 
field with periodic boundary conditions |36j : 



M 

H = Y. + ^I^I+i + ^i^l^Ui - 1)) - hSz, (1.1) 
fc=i 

where 

1 ^ 

^^ = 2^^^' [H,S,]=^. (1.2) 

k=l 

Here A is the anisotropy parameter, h an external classical magnetic field, and the length M 
of the chain is assumed even. The quantum space of states is = (S'^^Wfe, where Tik ~ is 
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called the local quantum space at site k. The operators c^'^'^ act as the corresponding Pauli 
matrices in the space 7ik and as the identity operator elsewhere. 

The description of the full dynamical properties of such a model amounts to the determination 
of its correlation functions and associated structure factors. This was first achieved in the case 
equivalent to free fermions, namely at A = for the XXZ chain. We stress that these free 
fermion results already required tremendous efforts to get satisfactory answers, in particular 
for the spin-spin correlation functions [73l[Mlll[2nilH[85l|72l[76l[7ai4l[7ll[2^ 
Going beyond this free fermion point turned out to be very involved. The first attempts to 
compute the correlation functions at arbitrary anisotropy A came from the algebraic Bethe 
ansatz [39] (see also [38] for the one-dimensional Bose gas). They led however to rather implicit 
representations. This was due to the complicated combinatorics involved in the description of 
the Bethe states, in the action of the local operators on them and in the computation of their 
scalar products. The notion of dual fields (which are quantum operators) was introduced in 
[58] to try to overcome such combinatorial difficulties. In this approach (see |59] and references 
therein), the correlation functions are obtained in terms of expectation values (with respect 
to the dual fields vacuum) of Predholm determinants depending on several such dual fields. 
Unfortunately, these auxiliary quantum fields can hardly be eliminated from the final answers, 
making the use of such results quite problematic. More explicit representations were obtained 
later on in terms of multiple integrals for the elementary blocks of the correlation functions; it 
uses a completely different approach, working directly with infinite chains (hence with several 
hypothesis), based on the representation theory of quantum affine algebras and their associated 
g- vertex operators [MlllHlllS]. More algebraic representations of these elementary blocks have 
been obtained recently along similar lines together with a link to fermionic operator expressions 
of the corresponding correlation functions [HI [131 [121 [15] . Extensions of this scheme to non zero 
temperature have been considered in [311 [32] . 

The resolution of the so-called quantum inverse scattering problem for this model (namely 
the reconstruction of the local spin operators in terms of the operator entries of the monodromy 
matrix) opened the way to the actual computation of the form factors (the matrix elements of 
the local spin operators in the eigenstates basis) and correlation functions [52], [53l [6^ in the 
framework of the algebraic Bethe ansatz. With the help of previous results (like determinant 
representations for the partition function with domain wall boundary conditions, and for the 
norms and scalar products of Bethe states [30l [571 EHl [371 EE]), this approach led to explicit 
determinant representations of the form factors of the Heisenberg spin chains in a magnetic 
field as well as to multiple integral expressions for the elementary blocks (or reduced density 
matrix) for their correlation functions in the thermodynamic limit. Remarkably, these integral 
representations coincide (in the zero magnetic field limit) with those obtained previously directly 
for the infinite chains using g-vertex operator methods and the representation theory of the 
quantum affine algebras [HJ [ISJ [12] • 

These advances immediately raised the next challenge, namely the question of obtaining 
manageable expressions for physical correlation functions such as the spin-spin correlation func- 
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tions. This is a central problem in the field of integrable models, on the one hand from the 
viewpoint of effective applications of these models to the realm of condensed matter physics 
(the associated dynamical structure factors are measurable quantities), and on the other hand 
concerning the more fundamental question of the long distance asymptotic behavior of the cor- 
relation functions. While the practical determination of the dynamical structure factors can 
effectively be addressed using the above form factors expressions in terms of determinants com- 
bined with powerful numerical techniques to sum up their corresponding series (see |20l I19|). 
answering the second question from purely analytical techniques looks like a fantastic theoretical 
problem. It amounts to understanding how the microscopic nearest neighbors interactions, for 
example in Heisenberg spin chains, integrate to produce an effective long distance behavior of 
the spin correlations. In particular, for models having no gap in the spectrum, the spin-spin 
correlation functions are believed to have, besides their possible trivial constant value, a power 
law decay with the distance. For generic integrable models, in particular not equivalent to free 
fermions, extracting such an asymptotic behavior "from first principles" has been a challenge for 
many years. 

In the case of fully interacting gapless models (like the XXZ Heisenberg spin-| chain for 
anisotropy — 1 < A < 1), the first predictions for this power law behavior came from a conjec- 
tured correspondence between such integrable lattice models with nearest neighbors interactions 
and a continuum theory having long-range interactions, the Luttinger model. Using this hypoth- 
esis, Luther and Peschel [6^ succeeded to predict the exact value for the XXZ critical exponents 
at zero magnetic field as functions of the anisotropy parameter A. The above correspondence 
was further enlighten by Haldane's development of the Luttinger liquid concept, leading in turn 
to predictions for the critical exponents of the XXZ model in a magnetic field [33], [Ml [35] . The 
Bethe ansatz techniques were used there to compute the parameters describing the low energy 
excitation spectrum of the XXZ model and to show that they indeed satisfy relations that are 
characteristic of the Luttinger liquid universality class. This conjecture was further studied in 
a series of papers [9j [TOl [40] where, in particular, the long-distance behavior of the longitudinal 
spin-spin correlation function was computed in perturbation to the second order in A around 
the free fermion point A = 0, with results in agreement with the above predictions. 

Another approach to this problem stemmed from the hypothesis that critical statistical 
systems with short range interactions should be described by a conformal field theory. In this 
framework, the mapping of a conformal field theory defined on the plane into one on a strip of 
finite width I led to the determination of the critical exponents in terms of the eigenvalues of the 
transfer matrix along the strip [III [U [ISj [8] . In particular, for large width I, the central charge 
and the scaling dimensions of scaling operators of the theory are given by the behavior of 
respectively the ground state and the different excited states energy levels. The possibility to 
compute the finite size corrections to the spectrum of integrable models in the framework of 
Bethe ansatz methods gave rise to a prediction for the corresponding critical exponents, in full 
agreement with the above bosonization approach [261 ESI [EHl [Ml [EH [SH [SSI [56l [25] . 

Further works in these directions (based however on several conjectures) sharpened this 
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picture by providing, in addition to the critical exponents, predictions (at zero magnetic field) 
for the amplitudes of the spin-spin correlation functions asymptotic behavior, which are not 
directly accessible from the bosonization or conformal field theory approaches [Ml [65l [66] , 



1.2 Overview of the method 

As already mentioned, the purpose of the present article is to describe an effective method 
allowing us to derive the long-distance asymptotic behavior of the correlation functions "from 
first principles". We will use the algebraic Bethe ansatz framework and our previous results 
|48l I45j giving an exact (master equation) formula for the spin-spin correlation functions of 
a finite XXZ chain. Starting from this formula, we will explain how to take the thermodynamic 
limit and extract the large distance asymptotic behavior of the correlation functions step by 
step from their lattice expressions. As the master equation representation also applies to other 
integrable models [IS], we will briefly describe how to implement the XXZ related analysis to 
the computation of the asymptotic behavior of the density-density correlation functions in the 
interacting Bose gas in one dimension. Let us now briefly explain the main features of the master 
equation representation (we refer the reader to the original articles [501 145] for more details). 

In the case of physical correlation functions, the dependence on the distance m appears in 
the form of the m*'^ power of the elementary shift operator. This operator in its turn, is some 
function of the transfer matrix of the integrable model at hand. The first step towards the master 
equation is to consider an integrable deformation (twist) of the transfer matrix (and hence of the 
shift operator) depending on some continuous parameter k. This parameter is chosen in such a 
way that one recovers the usual correlation function at say k = 1, and such that the spectrum of 
the twisted transfer matrix for k around some kq is simple and disjoint from the one for k = 1. 
These properties allow us to sum up exactly the form factor type expansion of the K-deformed 
correlation function with respect to the eigenstates of this «;-twisted transfer matrix (for any k in 
some neighborhood of kq) in terms of a single multiple contour integral that we call the master 
equation representation. The original correlation function one started with (corresponding to 
the point k = 1) is then reconstructed at the end of the computation from the knowledge of its 
values in an open neighborhood of kq. 



Originally, such a representation was first derived in [50j for the generating function of the 

longitudinal spin-spin correlation functions of the finite XXZ chain. This generating function 




is given by the ground state average value of the operator e^^™ [39], where (3 is an arbitrary 
complex number and 




(1.3) 



n=l 



SO that we have 
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(1.4) 



/3=0 
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where the symbol stands for the second lattice derivative : 

Dl f{m) = f{m + 1) + f{m - 1) - 2/(m). 
Setting K = , the operator e'^®'" can be rewritten as 

-''^'"=n(^+^<)' (1-5) 

and, thanks to the solution of the quantum inverse problem, it is shown to be equal to the m^^- 
power of the K-twisted transfer matrix multiplied by the m*'*-power of the inverse of the usual 
transfer matrix. Hence it is the simplest example for which the master equation representation 
can be obtained. 

Analogous representations were obtained for other correlation functions of the XXZ chain in 
|51| . Moreover, the time dependent case can even be considered along the same lines [49]. Later 
on, we also derived similar representations for correlation functions of a wide class of integrable 
systems (including continuous theories) possessing the i?-matrix of the six- vertex model [45j. 

One of the first problem to solve in this approach is that it is a priori not easy to take the 
thermodynamic limit directly in the master equation, for example for (e^^™). This is due to the 
fact that, in this limit, the number of integration variables involved into the master equation 
tends to infinity with the size of the lattice. However, starting from this representation, one can 
obtain various series expansions for the above master equation representation for the generating 
function. In these series representations the thermodynamic limit can be taken rather easily 
[501 is]. 

In this article we give a new expansion of the master equation, appropriate for the study of 
the long-distance asymptotic behavior of the generating function (e^^™). The key idea is to 
consider the master equation written as a multiple contour integral such that the only poles of 
the integrand inside the contour are at the Bethe parameters characterizing the ground state. 
It turns out that the algebraic structure of these poles is given by the square of a Cauchy type 
determinant, in complete analogy with the free fermion case. In fact this feature is quite general 
and will appear for other correlation functions and other models as well, although with various 
minor modifications. 

Then, the next step of the asymptotic analysis also takes its roots in the free fermion case. 
There, the obtained series is the expansion of a Fredholm determinant of an integral operator 
I + Vo (see ^M), 

oo 

(e^2-)^.o = E / detnVo{X„Xk) , (1.6) 

its kernel Vq depending on the distance m. It is known from the analysis of such kernel that 
the main contribution to the m — > oo asymptotic behavior is generated by the traces of powers 
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of the kernel Vq, tv{VQ). Independently of the power k involved, every trace behaves as 0{m). 
Therefore the natural idea would be to reorder the series p.6|l in such a way as to obtain an 
expansion with respect to tT{VQ). This can be easily done by presenting the determinants in 
(II. 6p as a sums over permutations provided the latter are ordered with respect to their cycles 
decomposition (see for example 

In the case of a general anisotropy parameter A, the series of multiple integrals for the 
generating function cannot be reduced to the form (|1.6p . Nevertheless, it can still be reordered in 
a way similar to the case of free fermions. Such a reordering leads to the appearance of multiple 
integrals of a special type that we call cycle integrals (see Section [3]). These cycle integrals 
play the role of analogs of tT:{VQ ). Each of them has a computable long-distance asymptotic 
behavior that we recently obtained using Riemann-Hilbert techniques applied to the Predholm 
determinant of an integral operator with generalized sine kernel L46J. Our strategy is then rather 
simple: using the asymptotic behavior for each cycle integral, we can sum up asymptotically 
the multiple series corresponding to the generating function {e^^"^). It so happens that such 
asymptotic series exponentiates in a natural way that mimic a Fredholm determinant expansion. 
However, in our case, all integrals are coupled in a non trivial way, and the correlation function 
itself is not a Predholm determinant (except in the free fermion point). We believe that this 
strategy can be applied to quite general cases although the details should be adjusted accordingly. 

As a result of this procedure, we are able to find the desired long-distance asymptotic behavior 
for (e^2™^. 

(e/32m^= GW(/3 + 27ri(T,m)[l + o(l)] , (1.7) 

(T=0,± 

where we have defined the function 

(/3, m) = C(/3) • e^^^^m^^ . (1.8) 

Here C{j3), D and Z are constants depending on the thermodynamic quantities that can be 
computed in terms of the low energy spectrum properties of the XXZ chain. This spectrum can 
be characterized by the solutions of the system of Bethe equations (|2.2|) [HITS]. In the thermody- 
namic limit, the ground state of the model appears to be the Dirac sea in the momentum space: 
the spectral parameters of the particles occupy the interval [—(7,(7] with a density p{X). The 
Fermi boundary q of the interval depends on the anisotropy parameter A and on the magnetic 
field h. Note that q remains finite for a non zero magnetic field, while it tends to infinity in the 
vanishing magnetic field limit. In the present article, q will be kept finite in all computations. 
However, the final expressions that we obtain are well defined in the zero magnetic field limit as 
well. In the Bethe ansatz framework, q and p{X) can be obtained in the following way. Let e(A) 
be the dressed energy, that is to say the solution of the integral equation 
q 

e(A) + ^jK{\-p) e{p) dp = h - 2p',{X) sin C, (1.9) 

-<? 
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where 

sin IC 

i^(A) = — — — — — , cosC = A, 0<C<vr, (1.10) 

smn(A + iQ sinn(A — iQ 

and "p^iX) is the bare momentum of the particle (see (12.6(1 ). Then q is such that e(g) = 0. 
Similarly, the ground state density p(A) [BOl E2] satisfies a linear integral equation 

P(A) + ^ 1 K(A - ^) p(^) = i-p'o(A), (1.11) 
-<? 

The dressed momentum is closely related to the density as 

KA) = 27r / p(^)d^. (1.12) 

JO 

The Fermi momentum, denoted as p^, is given by its value at the Fermi boundary, ■pp = p{q). 
The average density D is related to the ground state magnetization (cr^) and is expressed in 
terms of the Fermi momentum as 



{a^) = l-2D, D = I p{^,)d^, = ^. 



(1.13) 



Note that D = ^ in the case of a zero magnetic field hence implying a vanishing magnetization. 

Another constant appearing in the asymptotic behavior (|1.8|) is the quantity Z. It can be 
identified with the value of the dressed charge Z{X) at the Fermi boundary, Z = Z{±q). We 
recall that Z{X) satisfies an integral equation similar to (jl.lip 

<? 

Z{X) + ^ y K(A - /z) Z(/x) dfi = 1, (1.14) 

-<? 

and can be interpreted in the XXZ model as the intrinsic magnetic moment of the elementary 
excitations [59] . 

Finally, the constant C(/3) in (II. 8|) can be computed in terms of a ratio of four Fredholm 
determinants. The associated integral operators have compact support on contours surrounding 
the interval The explicit formulas are given in Section ISTTl At this stage, we would like 

to make two important remarks. 

The first one concerns the fact that all the Fredholm determinants involved in our compu- 
tations, for example those appearing in the constant C(/3), are associated to bounded integral 
operators acting on some compact contours surrounding the finite interval [—q,q]- Due to the 
finiteness of g, all these determinants take finite values, which is a direct consequence of working 
at a non zero value of the magnetic field. When g — > oo the contours are no more compact 
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and the corresponding Fredholm determinants diverge. It can nevertheless be checked that, 
in the zero magnetic field limit, the overall constant C(/5) appearing in the final answer (jl.Sp 
remains finite although each of its individual parts (i.e. the four above mentioned Fredholm 
determinants) diverges in this limit. 

The second remark is the following: since the operator e'^^" (II. 5p is a polynomial in k = , 
the generating function is a 27rz-periodic function of f3. Of course, the leading term of the long- 
distance asymptotic behavior is itself not necessarily a periodical function of /?. Nonetheless, our 
analysis shows that this a priori broken periodicity is at least partly restored by the first oscillating 
correction to the leading term, as it can be seen from the equation (II. 7p . More precisely, the 
result of the direct computation of the first oscillating correction to the asymptotics (|1.7p can 
be reproduced by a simple shift /3 — > /3 ± Ini from the non-oscillating term. This property, that 
we will prove here, has been recently used in p[6j to get predictions for the asymptotic behavior 
of (e^®™) at zero magnetic field. 

The long-distance asymptotic behavior of the longitudinal spin-spin correlation function then 
can be extracted from the one of (e^^™) using (|1.4|) . At the first leading orders, it reads, 

('^l(^m+l>,eading = (^^ ) " + 2 1 | ' ^2^2 , (1-15) 

where is related to the properly normalized form factor of the operator between the 
ground state and an excited state containing one particle and one hole sited at the two different 
boundaries of the Fermi sphere, namely at q and —q; hence, it corresponds to an umklapp process 
[27l [351 [34] . The factor 2 in front of |i^crP just corresponds to the two possible states of this 
kind. This result agrees with bosonization and conformal field theory analysis. Note however 
that, in the thermodynamic limit, all form factors scale to zero as some power of the size of 
the lattice. The precise coefficient and proper scaling behavior of the form factor which results 
from our computations is therefore rather non-trivial. We have nevertheless checked that, at 
zero magnetic field, this constant goes to a finite value, and that it agrees to second order in 
A around the point A = with the predictions given in [Ml [65l [66] . 



This article is organized as follows. 

In Section [2l we explain how to obtain, in the case of the XXZ spin chain, a new series 
expansion for the master equation representation of (e^^™), that is suitable first for taking 
the thermodynamic limit and then for extracting the asymptotic behavior of this correlation 
function. 

This expansion is written in Section [3] in terms of cycle integrals. Then, we recall the results 
of [l6] and apply them to get an asymptotic expansion of the cycle integrals in a form adapted 
for their asymptotic summation. 

The asymptotic summation of these cycle integrals is described in Section HI It is divided 
into several steps according to the different nature of the terms to sum up. In particular one of 
the steps uses a generalization of the Lagrange series (the details are given in an appendix). 
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This procedure enables us, in Section [5l to describe the asymptotic behavior of the cor- 
relation function (e^^™) at large distances m. We show the existence of both oscillating and 
non-oscillating terms that we compute at the leading orders. The asymptotic behavior of the lon- 
gitudinal spin-spin correlation function follows and is also given in the same section. It confirms 
at the leading order the predictions from bosonization and conformal field theory. 

In Section [6l the whole procedure is applied to another model, the one-dimensional interact- 
ing Bose gas, for which we compute the asymptotic behavior of the density-density correlation 
function. 

Finally, various technicalities are gathered in the appendices. We would like in particular to 
draw the reader's attention to Appendix [C] which contains rather essential ingredients for the 
summation described in Section [H 



2 Master equation and its thermodynamic limit 



As already mentionned in the Introduction, the main object of study of this article (at least up 
to Section [5]) is the ground state expectation value of the operator e^^"^ (|1.3p for the XXZ model 
(II. 1|) in an external magnetic field. In this section, we recall the master equation representation 
obtained in [501115] for this operator in the finite chain, and show how such a master equation can 
be expanded into some series suitable both for the thermodynamic limit and for the asymptotic 
analysis performed in the following sections. 



2.1 Master equation for (e'^^'") 

The master equation gives a multiple integral representation for the expectation value of the 
operator e^®" with respect to an A^-particles Bethe eigenstate |'0({'^})) of the Hamiltonian. It 
reads. 



(-1) 

A^! 



(V'({A})|V'({A})) 



r({A}) 



n 



dzj l{zj)d{zj) 
l{Xj)d{Xj) 



det^r?«({z},{A}|{z}) 



N 



(2.1) 



The various quantities entering this representation are all rather universal objects in the context 
of the Bethe ansatz. Let us explain their meaning. 

The parameters Ai,...,AAr describe the specific eigenstate \tp{{X})) in which the average 
value is computed. Prom now on, this state will be the ground state of the XXZ chain p.ip 
in a given magnetic field h and for some anisotropy A = cosC, < C < tt. In this case, these 
parameters are real numbers which satisfy the system of Bethe equations 



3;(A,|{A}) =0, j = l,...,N, 



(2.2) 
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with 

N N 



3^(/x|{A}) = a(/x) J] sinh(Afc - /i - iC) + d{fi) J] smh(Afe - /x + . (2.3) 
fc=i fc=i 

Here a(/i) = a^{fj,) and = d^{fj,) are the vacuum eigenvalues of the transfer matrix, with 

ao(/u) = sinh(^ - (io(/i) = sinh(/z + ^). (2.4) 

The whole explicit dependence of (|2.ip on the distance m is contained in the function l{z), 
which can be expressed in terms of oq and do as 



.MfJ-)/ ysinh(^+§ 

where 

(sinh(^ — /i) 
sinh(f + /i) 



Po(/u) = -ilog ( . ^ , ) , Po(0) = 0, (2.6) 



plays the role of the bare momentum of the pseudo particles (see also p.lOj) . p.lljl ). 

Besides the function 3^(Aj|{A}), the master equation (12. ip contains also a function y^izjllz}), 
K = , which naturally appears in the definition of the twisted Bethe equations associated to 
the spectrum of the K-twisted monodromy matrix (see [50] )• It is given as 

N N 
3^k(/"|{^}) = a{fi) sinh(2;fc + Kd{fi) JJ sinh(2;fc - fi + iC) . (2.7) 

fc=i fe=i 

One more object entering (|2.ip is the determinant of the matrix 0^, with matrix elements 

N 



i^f<.)jk ({^}' {-^}l{^}) = t{zk, Xj) JJ sinh(za - Xj - iC) 

a=l 

N 

- Kd{Xj)t{Xj,Zk) sinh(2;a - Aj + iC) , (2.8) 



a=l 

where 

— isinC 
smh(z — A) smh(z — a — iQ 

The determinant of this matrix describes the scalar product between any eigenstate of the twisted 
transfer-matrix and some arbitrary state (see [501 145] ). 

It remains to describe the integration contour r({A}) appearing in equation (12. ip . For every 
variable zj the integral is taken over a closed contour surrounding the points Ai, . . . , Aat and such 
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that any other singularities of the integrand (i.e. the roots of the system 3^K(zj|{2;}) = 0) are 
located outside this contour. Since here the parameters Ai, . . . , Xjy decribe the ground state of 
the XXZ chain in a non-zero magnetic field h, they are real numbers and, in the thermodynamic 
limit, they fill the interval [—q, q] (with q finite, see (11.91) ) with a density p{\) satisfying the linear 
integral equation (II. lip . Therefore, the integration contour will be chosen such as to surround 
the interval [—q, q]. 

It is worth mentioning that the integral representation (|2.1I) only holds for |fi;| small enough. 
This restriction does not cause any problems: the expectation value of e^^" is a polynomial 
in k; thus it is enough to know this polynomial in a vicinity of k = and then consider its 
analytic continuation. One should however remember that it is not always possible to analytically 
continue directly the integrand of (|2.ip . For instance, at k ^ 1, one of the solutions of the 
system yKizj\{z}) = goes to {A}, and an integration contour r({A}) satisfying the above 
desired properties does not exist anymore. 



2.2 Transformation of the determinants 

The integrand of (|2.ip has poles at Zk = Xj and, formally, the multiple contour integral in (|2.ip 
can be computed by residues at these points. These poles are contained in the determinants 
detr^K- More precisely, each of these determinants has simple poles at Zk = Aj, which means 
that the integrand has double poles at these points. 

In order to be able to analyze the contribution of these poles, it is convenient to extract 
them explicitly from detf^K. This extraction can be done in several ways. The method we 
present here consists in extracting a Cauchy determinant from each det using the fact that 
the parameters A are solution of the Bethe equations. It allows us in particular to show that 
det is proportional to k — 1. 



Proposition 2.1. If the parameters Ai 
detN^Ki{z},{\}\{z}) = detAT 



, Xn satisfy the system of Bethe equations, then 

N 



1 



1 



vr\0) - kVi\9) 



sinh(Aj — Zk) 

N 

sinh(2;a 

a,b=l 



n \ «(^^) 



V_(A 



1 



Xb - iC) ■ detN (Sjk + U^ki<^)) . (2-10) 



or 



detN^ni{z},{X}\{z}) = detAT 



1 



sinh(zfc — Xj) 



N 



N 

n \ ^(^^■ 



V-{e)-KV+{e) 



II sinh(Aa -Zb- iC) ■ detN (6,k + f^jfcV)) • (2.11 



a,fe=l 
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In these expressions, 6 is an arbitrary complex number, and 

^ sinh(/i — Aa ± i() 



{A} _ 

{•2} J smh(/i - Za± iC) ' 



n 



(2.12) 



The entries of the matrices U^^\6) and U^^\9) are given by 



N 

n sinh(za 

a=l 


-A,) 




-Afc) 


-K^{9- 


Xk) 


N 

n sinh(Aa 


-A,) 


y+ 


'(A,)- 




) 


a=l 












N 

n sinh(zfc 

a=l 






- ^;fc) - 


- Kf,{zj - 




N 

n sinh(zfc 


- Za) 






KV+{Zk) 




a=l 













(2.13) 



(2.14) 



with 



(2.15) 



K^{X) = coth(A + iC) - Kcoth(A - <) . 

The proof of this proposition is given in Appendix El 

Remark 2.1. One can observe that detri^ is proportional to 1 — k. However, the second line of 
(I2.10p or of (|2.11|) does not necessarily vanish at k = 1 (/? = 0). Indeed, if Zj = Xj (modulo 
the permutation group) for all j = 1, . . . , then V±{6) = 1 and the factor 1 — k cancels. Such 
situation can arise in the process of evaluating the multiple integral (|2.ip . 

In the representations (I2.10p . (12. lip , the poles at Zf^ = Xj are gathered explicitly in the 
two Cauchy determinants. The entries of the matrices U^'^'^^ contain singularities at Xj = Afc 
(respectively Zj = Zk), but the corresponding determinants are obviously not singular, since 
the original det fi^ vanishes at Xj = Xk or Zj = z^. To make this fact explicit, we present 
the determinants of the finite size matrices 5jk + Uj^'^^ as Predholm determinants of integral 
operators acting on a contour T surrounding the points {A} and {z} (see Fig. 1): 



det 



N 



S3k + ui','\o) 



det 



W, w') 



where 



N 

n 



sinh(tt; — Za) Kt^{w — w') — Kk{6 — w') 



^}-^ s\n\i{w - Xa) V^^{w) - nVl^{w) 



and 



N 



sinh(ti;' — Aq) Ki^{w 



w 



sinhfw' — Za) V-(w') — kV+(w') 

a=l 



(2.16) 



(2.17) 



(2.18) 
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r 




—q Ai A2 


X .1 


r • x X 









Figure 1: The contour V . 



The equivalence of the two representations (|2.16|) is proven in Appendix IA.2[ 

Using these results, we can rewrite the master equation (12.11) in the following form: 



N 



f n 



r({A}) 



27ri 



'^^ y ^) d(Zj) l}^ smh(zj-2a-iC) 



1 



Ai 



det 



1 



N 



sinh(zfc — Aj) 



(2.19) 



where 



Wn 



Ai, 



^ / {A} \ _ -TT sinh(2a - Afc - iQ sinh(Afe - Zg - iC) 
ZN j \ {z} j ^^^^ sinh(za - Zb - iQ sinh(Aa - A^ - iQ 



{v+\9i) - KVi\ei)){v.{e2) - KV+{e2)) 



(2.20) 



and 



e 



jk 



N \ ^1 

KAj) W sinh(Aa - \j - K) 

a=l J 



dy{Xj\{X]) 
dXk 



(2.21) 



It is important to note that the function W (|2.20p is a symmetric function of the variables {z} 
and also of the variables {A}. Moreover, it possesses the recursive reduction property 



Wn 





• , AAr\ 




ZN J 



Wn- 



Ai 



A; 



Z]V=Ajv 



, . . . , ..N-1 
Zl, • • • , ZN-1 



(2.22) 
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The entries of the matrix contain the Lieb kernel K(\) defined in (ll.lOp . 



Qjk = Sjk I log 



(2.23) 



2.3 Expansion of the master equation 



We now express the master equation as series appropriate for taking the thermodynamic limit. 
Replacing one of the Cauchy determinant by the product of its diagonal elements we re-write 
equation (|2.19|1 as 



det/vT 



N 



r({A}) 

N 



{A}^ 



1 



k=i 



sinh(zfc - Afc) 



det 



N 



where 



Vat /i 



kV+\ 



VI' I ^ 



hNizk\{z}\{X}) - 1 
sinh(zfc — Aj) 



{A}^ 



+ 



Wn 



1 



{A} 
{4 



sinh(2;fc — Aj 



, (2.24) 



(2.25) 



and 



K-V- { Zk 



hNizk\{z}\{\}) 



{A}^ 



' ^ ) d{zk) 11 sinh(zj.-Za-0 
a=l 



(2.26) 



Note that the function V defined by (|2.25|) satisfy similar properties to W: it is also a symmetric 
function of the variables {A} and of the variables {z} separately, and possesses the same reduction 
property: 



Vat ^JL 



Ai, . . . , \n 

Zi, ZM ^ 



Vn-1 fJ- 



Ai, . . . , Aat-i 

Zi, . . . , Zjy^i 



(2.27) 



•2jv=Ajv 

As for /iTv (|2.26p . it possesses also interesting reduction properties. Indeed, the Bethe equations 
(12.21) for {A} ensure that, for any k and s in 1, . . . , A^, 



hNizk\{z}\{\})\ 



Zk=\s ' 



(2.28) 



and 



d 

dzk 



hN{zk\{z}\{\}) 



V 



N-1 1 



A. 



{A}\A, 

{z} \ Zk 



1 / q(a^) . 



N \ 

Y,K{Xs-Xa)]- (2.29) 

a=l ) 
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It follows in particular from (|2.28|) that the first part of the determinant in (|2.24|) is holomorphic 
at Zf^ = A J . 

Let us now expand this determinant via the Laplace formula. We obtain 



1 



detpj B 




r({A}) 



N 



^ ^ (-l)[^(-)]+[^W]detfc 



Jm{po{zj)-po{Xj)) 

sinh(zj — Xj) 

hNizk\{z}\{X}) 




Wn 




n=0 {z}={z}.yU{z}^ 
{X}=_{\}cU{X}^ 
#7=#a=" 



G7 



sinh(2fc — Xj) 



detfcg^ 



sinh(z,fc — Xj) 



(2.30) 

The above sum runs through all possible partitions a L) a and 7 U 7 of 1, . . . , A'" such that 
#7 = #a = n. The sets of parameters {A} and {z} are divided accordingly into disjoint subsets 
{A} = {A}q,U{A}q,, and similarly {z} = {z}^U{z}^. Moreover, we specify that the parameters in 
each subset are ordered in the canonical way, namely {A}q, = {Xai, Xa2, ■ ■ ■} with ai < 02 < • • • 
etc. Finally, [-P(a)] (respectively [-P(7)]) stands for the signature the permutation P such that 
P(q, a) = l,...,N (respectively ^(7, 7) = 1, . . . , A^). 



Lemma 2.1. In the multiple integral 112. 30\) . 

1. the only non-zero contributions come from partitions such that 7 = a; 

2. the first determinant in 112. 30]} contributes only through the product of its diagonal elements. 

Proof — Let us consider a partition such that 7 7^ a. Then there exists £ G 7 such that i ^ a. 
Let us consider the integrand of (|2.30p as a function of zi. At first sight it has a simple pole at 
Zi = Xi and also simple poles at the points Zi = Xj for all j € a. Yet, due to the property (|2.28|) . 
all the corresponding residues are equal to zero. Hence, the integral over zi vanishes. Thus we 
conclude that only partitions such that ■y = a yield non- vanishing contributions to the integral. 

It remains to prove the second part of the lemma. Let us consider an off-diagonal contribution 
coming from the first determinant. Hence there is a variable Z£ G {z}'y such that the integrand, 
considered as a function of this Z£, has at most simple poles. However, due to the presence of 
the factor hj\f{zi\{z}\{X}) — 1, their residues vanish. □ 

Thus, 



r({A}) \ ^ V ; 




N 



S S S ^i-M-. - A.) n 

#a=n 



hN{zk\{z}\{X})-l 

smh^{zk - Xk) 



det 



'j,kS:a 



sinh(zfc — Xj] 



(2.31) 
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The integrals over {z}^ can now be easily computed using (I2.29p . We obtain 



/3S„ 



detAT©^ ^ / I 27rz \ {zja / I 

#a=„r({A}) 

( j'l: ) ■ n (2":Mp(A,,) . n -i;^ ■ <ie...« -i^. (2^32) 



where 

27riMp(A, ) = log' ^ - i V /^(A, - A„) . (2.33) 



The notation p(A) is motivated by the fact that p{X) goes to the ground state density p{X) (jl.lip 
in the thermodynamic limit. 

It remains to replace the summation over partitions by a sum over individual A's. Let us 
denote by A the original set of spectral parameters {Ai,...,AAr} describing the ground state. 
We then have 



/3Q„ 



nf=i (2vriMp(A,)) 1 



n=0 Ai,...,A„gA j=l ^ '^^ 

-|--^ sinh(2;j — Xj) " \^ {2;} ^ " sinh(zfc — Aj 



n d^^- fr e-(^ofe)-^o(^.)) / {A} I ^ 



X 

r({A}) 

with 

' {A} ^ _Tr. i {A} ^ jj„ , {A} 



2.4 Thermodynamic limit 

Recall that the thermodynamic limit corresponds to A^, M — > 00 in such a way that the average 
density D = N/M remains fixed. In this limit, discrete sums over the ground state Bethe 
parameters Afe turn into integrals over [— g] weighted by a density function /o(A), 

^ E /(A,) - / /(A)p(A)dA. (2.36) 

It is easy to compute the limit of the pre-factor in (|2.34p . Using the notation (|2.33p . we can 
express the determinant of the matrix G (I2.23P as 

N 



detAT e = {2mMp{\^)) ■ deiN 



, K{Xj - Afc) 
^''^ 27TMp{X,) 



(2.37) 
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The equation (|2.33|) for p{X) can be written in the form 



1 ^ 



(2.38) 



a=l 



Replacing the sum by the integral and using the integral equation for the ground state density 
(jl.lip . we find that in the thermodynamic limit /5(A) — > /9(A). Therefore the determinant in the 
r.h.s. of (|2.37|) goes to the Fredholm determinant of the integral operator I + acting on 
[-q,q], namely, 



N 



detTv e • JJ {27riMp{Xj)y^ det 



I+—K 
27r 



(2.39) 



where the integral operator acts on the interval [—q,q]- 

Hence, in the thermodynamic limit, the generating function (e^^"") is given by the following 
series of multiple integrals: 



1 °° 1 

det (/ + ^K) ^ n\ 

\ l-K I n=0 



r(hg,g]) 
{A} \ 

{A ] 



i=i 



2m 



" eim{po(zj)~po(\j)) 



sinh(2;fc — Aj) 



. (2.40) 



The question of convergence of this series is quite complicated. To understand this topic a 
little better, let us consider the free fermions point A = (C = f )• 

In this point the terms of the series can be significantly simplified. Indeed, it is obvious from 
(I2.12P that Vn = K — 1 for C = f • It is less trivial, but true (see Appendix IbI). that Wn = 1. 
Therefore Tn = {i^ — 1)", and we have 



E 

n=0 



(k- 1)" 



li Otti / 11 '}.7ri 



r([-9,q]) 



-1 ^ = 1 



" gjm{po(zj)-po{Aj)) 



1 



sinh(2;fc — Aj 



(2.41) 



The integrals over are now factorized and can be taken explicitly. We obtain 



(k - 1)" 



n=0 



n! 



dAj • det„ 



sin(f(po(A,)-po(Afe))) 
vr sinh(Aj — \k) 



(2.42) 
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The series (|2.42p is an expansion of the Fredholm determinant of the integral operator / + Vq 
with kernel 

( ,^ sin{f [po(A)-po(^)]} 

VbA,^ = K-l) — T ^ . 2.43 

7rsinh(A — ^) 

The general theory of Fredholm determinants ensures that the series (|2.42l) is absolutely con- 
vergent and define an entire function of n. It means that the expansion coefficients decay faster 
than exponentially. 

In the case of general A, the series (|2.40|) cannot be reduced to a simple form as in (|2.42|) . As 
mentioned above, the analysis of its convergence is therefore rather complicated. Nevertheless, 
taking into account that the XXZ chain with a general anisotropy parameter can be considered 
as a smooth deformation of the A = case, we shall assume in the following that the expansion 
(12.40p also is absolutely convergent. 

As a conclusion to this section, we would like to stress the role of the non-zero magnetic field. 
As we already mentioned, the support of the ground state density at non-zero magnetic field is 
a finite interval [— ^, Therefore, all the multiple integrals in the series (I2.40p are convergent. 
The pre-factor det[/ -|- -^K]'^ is also finite. However, at zero magnetic field, q ^ oo, and it is 
easy to see that each separate integral in the series (|2.40|) becomes divergent. This divergence 
should be of course compensated by the divergence of the Fredholm determinant det[/ -|- -^K]. 
The latter becomes ill-defined at g — > oo, because the kernel of the integral operator depends on 
the difference of the parameters, i.e. K = K[X — It is not clear, however, how one can extract 
the finite part from the representation (|2.40|) directly at /i = 0. Therefore one can say that the 
non-zero magnetic field plays the role of a regularization in the framework of our approach. One 
can check however that, in the end of the computation, the obtained results have smooth limits 
at /i — > 0, although the corresponding proof is highly non-trivial. 



3 Cycle integrals 

Even in the simplest case of free fermions where the multiple integrals are decoupled, the cal- 
culation of the large m behavior of (e^^™) remains non-trivial. Indeed, one has to analyze the 
large m behavior of the Fredholm determinant of the integral operator I + Vq whose kernel is 
given in (|2.43l) . In the general case, the situation is even worse as the multiple integrals are 
highly coupled and hence the corresponding series cannot be identified with the Fredholm series 
for the determinant of some integral operator. 

However, the series (|2.4Q|) can be rearranged in a specific way, making possible the asymptotic 
analysis of its separate terms. The key observation towards such a rearrangement comes from 
the free fermion point. There the leading asymptotic effect is produced by traces of powers of 
Vo : tr(V[f) behaves as 0(m) when m — > oo for s = 1, 2, . . . . Analogs of such traces of powers of 
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Vq will be given, in the general case, by the following type of integrals: 

' ^^'^ ~ J (2^ J (2^ I {z} j JJ- sinh(z, -A,)sinh(z, -A,+i) ' ^ 

r([~g,g]) -q ^ ^ 

in which the 2s integrals are coupled via some holomorphic function Qs of 2s variables Ai, . . . , A^, 
zi, . . . , Zs, symmetric separately in the s variables A and in the s variables z. Note that, in (I3.1|l . 
we identify A^+i with Ai. 

In this section, we explain how the representation (|2.40p can be decomposed into such inte- 
grals, and how the asymptotic behavior of the latter can be analyzed in the large m limit. 

3.1 The cycle expansion 

The decomposition of the series (|2.40p into multiple integrals of the form (13. ip . which plays 
a key role in its asymptotic analysis, comes directly from the cycle expansion of the Cauchy 
determinant. In this article we use the following definition for a cycle: 

Definition 3.1. Let A = (aij) be an arbitrary n x n matrix. A cycle of length £ is any product 
of entries of A of the form 



With this terminology, integrals of the form (13.11) will be called cycle integrals. 

It is well known (see e.g. jTl]) that the determinant of any matrix A can be presented as a 
sum of products of cycles of different lengths. We call such a representation of a determinant 
its cycle expansion. In our case, the determinant we need to decompose into cycles (namely, the 
Cauchy determinant in (|2.40|) ) is part, together with some symmetric function, of the integrand 
of a n-fold integral. For such kind of representations we have the following result: 

Proposition 3.1. LefWx and'^z be two curves inC, g a continuous function in {z, A) G x'^A; 
and Qn a continuous function in n variables A and n variables z on 'lo^ x "^J*. Assume more- 
over that, for any permutation a G S„, Qn is invariant under replacements of pairs {Xi,Zi) i-^ 
{K{i),Za{i))- Then, 

^Jd-Xj d". Qn ( ] det„ [g{z,,X,)] 

= C{nm) j d"A / d^z Qn ( J^J ) n 9{ztAt+i), (3.2) 

k = l 
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where 



c{nm) = n 



s=l 



V. 



s+l \ I 



(3.3) 



Here J^£j denotes the following set of triplets, associated to the configuration given by £i, . . . ,£n, 
J{e} = {{s,P,j) ■.l<s<n,l<p<£s,l<j<s}, (3.4) 
and, for t = {s,p,j), t+1 denotes the triplet {s,p,j + 1), with {s,p, s + 1) = {s,p, 1). 

Remark 3.1. The notations introduced above should be understood as follows: s labels the 

length of a cycle, ig stands for the number of cycles of length s and j marks the position of 
the variable in the corresponding cycle. Hence, for t = {s.p. j), fit = l-i's,p,j stands for the j^^ 
variable of the p^^ cycle of length s. In this notation we assume that iJ,s,p,s+i = l^s,p,i- 



Proof — Let ^ be a n x n matrix. Then, 



1 5" 

detnb] = - — detn[I + 75] 



1 & 



1 5" 
n! dj^ 

1 a" 



JJexp 



s=l 



-y=0 IT-- 



-exp(trlog[/ + 75]) 



7=0 



' ei,...,in=o S=l 

n n 



7=0 



7=0 



(3.5) 



s=l 



Substituting here gij = g{zi,Xj) and writing explicitly all [tr((?'')]^'', we obtain 



n is 



det„[5(^,,A,-)] = E^HW)nn E n^K^.r^toi))' M 

S=l p=l i^s,p,l) v,j(s,p,s)=l J=l 



^lv--)^n — 

Sfc£fc=n 



It remains to observe that, in the product of the functions 9{zi(^s ^ K^^^ p j+i)) ^ H^^PJ) should 
be different, otherwise the corresponding term does not contribute to the determinant (all other 
contributions eventually cancel out). Hence, the multiple sum over i(s,p,j) is in fact the sum over 
permutations: 



detn[g{Zi,Xj)] = ^ C{n\{i}) Yl 9{Za(t),Xa(t)+l)- 



(3.7) 



T:k£k=n 



Due to the symmetry properties of the function Qn, all the terms of the sum over permutations 
give the same contribution to the integral. n 
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Let us apply this result to the Cauchy determinant under the multiple integrals in the series 
(ITiO]) for (e^2m^ Denoting 



we obtain 



det [/ + i^K] 



(3i 



q 



n=0 £i,.../„=0 

SA;4=n 



d"A 

(27ri)" y (27ri)" 



*^ n. 



n 



{A} \ 

pim(po(2t)-po(At)) 



sinh {zt - \t) sinh (z^ - Aj+i) 



(3.9) 



It is actually convenient to remove the constraint Tiki^ = n on the summation variables tk in 
the sum (|3.9p . This can be done by introducing, in the n**^ term of the series, an n*^-derivative 
over some auxiliary parameter 7, similarly as in the proof of Proposition 13.11 



E 

n=0 



1 a" 

n! dY' 



00 "1 

E Ht^ 



,...,^„=0 A:=l 



nS/ n 



27ri 



X 



k{0 



{A} 

{4 



n 



r([-g,g])*e-^{0 - 
pim(po{2t)-po(At)) 



27ri 



sinh {zt - \t) sinh {zt - Xt+i] 



(3.10) 



7=0 



where \ J{i}\ = Efc^fc denotes the cardinality of the set J{£} (this cardinality is no longer equal to 
n since we have removed the constraint). 

The v}^ term in the series for (e'^®'")) is thus expressed as some sum of multiple cycle integrals 
of the form (13. Note however that these cycle integrals are coupled through the function 
which involves the total number 2|J|£j| of variables, and not only the 2s variables associated to 
a given cycle. In order to formalize this fact, we need to introduce some notations. 

Let us first recall the standard order on N^: 

{h,t2-,h) < {ui,U2,U3) <^ {ti < ui or {ti = ui and {^2 < U2 or {t2 = U2 and ^3 < us}}}}, 
and let it) be the identity operator on functions of one variable A and one variable z: 



id 



A 

z 



A 

z 



(3.11) 



For some cycle configuration {i}, we define I(s,p) as the operator acting on functions of 
2| J{£}.| variables Xt and zt, t € J{e}^ as the cycle integral Ig (13.11) over the variables Xg^pj, Zg^pj, 



24 



1 < J < i.e. 

= (8) ^ (8) ® ^ • (3.12) 

t<(s,m,l) variables t>(s,m,s) variables 

Note that the action on G\j^^y\ of a single 1(s,p) (associated to a given cycle (s,p)) still produces 
a function (with respect to the other variables). However, the action on (denoted here 

with the symbol *) of the product of such operators I(s,p) with respect to all cycles, 

fin/-) *^'--finn ) ( {t:i ) 

" -il- * g«m(P()(2s,P,j)-Po(As,p,j)) 

^ s=i p=i J=i '^^'^ ^^''P'j ~ ^'^'P'j) ^'""^ (^^'PJ ~ ^^-P'i+i) ' ^^^^^ 

produces a number, since all variables A and z are integrated. 

This notation enables us to recast the generating function in a quite compact form: 



n=0 



(3.14) 



7=0 



where 

~k / 

ti,...,e„=ok=i ^ \s=ip=i 



« _n I — 1 \ / \ „_i „_i / 



Therefore, the whole problem of the asymptotic analysis of the series (ITinll for (e^2™^ 
reduces to the derivation of an asymptotic formula for Gi...„(7). The advantage of such a 
representation is the following: in spite of the fact that all cycle integrals are coupled with each 
others through the function the asymptotic behavior for large m of every set of cycle 

integrals can nevertheless be computed separately, independently of the other integrals; then 
the asymptotic formula for multiple cycle integrals can be obtained by applying consecutively 
the results of the analysis of each pure cycle integral 1^(s,p) to the corresponding group of cycle 
variables {^s,p,j}i<j<s and {-Zs,p,j}i<j<si while considering the remaining variables fixed. Passing 
from the exact formula to the asymptotic one produces a very powerful effect: due to peculiar 
properties of the function the cycles become "quasi-decoupled" in the m +oo limit. 

Such a quasi-decoupling is enough to perform the sum over all the possible numbers £k of cycles 
of length k. 

The summation procedure of the asymptotics of multi-cycle integrals will be presented in 
Sectiondl However, first we need to remind the results of [46] concerning the asymptotic analysis 
of pure cycle integrals. 
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3.2 Asymptotic analysis of cycle integrals 

The asymptotic study of cycle integrals of the form (|3.ip was performed in 06] from the one of 
the Fredholm determinant of a generalized sine kernel. It was shown there that such integrals 
admit, up to any arbitrary order no, an asymptotic expansion of the form 

no -, 
Ul 

n=l 

+ E e-[-(^)--(-)l ± ^jt''m + 0('^). (3.16) 

rGZ* n=2\r\ ^ ^ 

The leading behavior (up to o(l) corrections) is given by /i'^'* which contains a linear term 
in m, a term proportional to logm, and a constant term. Its explicit form, which was computed 
in [l6], will be specified later. This leading term admits two kinds of corrections: oscillating 
(containing an oscillating factor e*'''"[P''^'^^~^°*^~'^^]), and non-oscillating oneJl]. Note that the 
leading term itself is non- oscillating. 

In order to derive the asymptotic behavior of (e^^™), one will have to sum up (|3.16p according 
to (|3.15ll . If we suppose that subleading corrections in (I3.16P remain subleading through this 
process of summation (which is not completely true, as it will be explained later), then it is in 
principle enough, so as to obtain the leading asymptotic behavior of (e'^^"), to consider only 
the leading term li^\Gs] in (|3.16|) . Recall however that our main goal is not (e^®™) itself, but 
the spin-spin correlation function (cf obtained from (e^^™) by second derivative with 

respect to (3 and by second lattice derivative (see (|1.4|) ). This process of taking the second lattice 
derivative will decrease the order of the non-oscillating terms, but not of the oscillating ones, as 
we may merely differentiate the exponential. Therefore, to be sure to get the leading behavior 
of (cf cT^_|_i), we will have to consider not only the leading (non-oscillating) term but 
also the leading oscillating correction, given by 

OslQs] = Otm + OiiGs], with of = -Le±^--lPoii)-Poi-i)]ii2;±i), (3.17) 

It will indeed appear that the non-oscillating term gives rise to an behavior for (erf a^_^_l), 
whereas the oscillating one produces an decrease, 6 being lower or higher than 2 according 
to the value of A. 

To specify the functional action of and Og on the function Gs and to apply it to our 
particular case, let us first briefly remind the main idea of the analysis performed in [46] . 

^li"'^\Qs] and /i"' [Ss] still depend on m, but simply as a polynomial (of degree at most s) of logm. 
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3.2.1 From generalized sine kernel to cycle integrals 

Suppose that the function Gs appearing in (13. 1|) is realized as a product of functions of one 
variable, i.e. 



{A} 
{4 



n^(A.)e^(^'), (3.18) 



1=1 



with g and f holomorphic in some vicinity of [—q, q]. Then the cycle integral ()3.ip can be 

explicitly computed in terms of the s^^ 7-derivative of the Fredholm determinant of the operator 
/ + on [—q, g], with kernel 

V{X,,) = F(A) ""^^ N(A) -PoM] - iiaW-am ^ F(A) = e^W,,(A). (3.19) 

7rsinh(A — fi) 

Indeed, from the series representation of log det + one has, 
s-l ^ s 

^~-^\., Oflogdet [I + ^V] = [d'XY[V{\k,Xk+i), with A,+i = Ai, (3.20) 
(s-ij! 7=0 ^J-^ 

which is exactly the cycle integral Is [Qs] once the residues in z have all been computed. 

In order to obtain the asymptotic expansion for cycle integrals involving more general holo- 
morphic functions Qs that are symmetric with respect to the variables Ai, . . . , A^ and zi, . . . ,Zs 
separately, one uses the fact that such functions admit an expansion of the form 



{A} 



Y,IlM^j)Mzj), (3.21) 

k=ij=i 



where (pkiX) and (pkiz) are holomorphic in a neighborhood of the interval [—q,q]. Setting then 
Fk{X) = v^fc(A)0fc(A) e^'^W = MX), (3.22) 
and denoting the corresponding kernel (I3.19P as 14, we obtain 

[^^] = E JTZrui ^7 log det [/ + jVk] . (3.23) 
fc=i ^ 

Hence, the asymptotics of cycle integrals for such Qs can be inferred from those of the Fredholm 
determinant of / + 'yV, V being a generalized sine kernel of the form (I3.19P (see |!46j for a more 
rigorous proof of all this procedure). 

The large m asymptotic behavior of log det [/ + 7!/] was obtained in [IH] by Riemann-Hilbert 
approach. There it was proved that, in the m +00 limit, 

logdet[/ + 7y] = >Vo(m) +0(1), (3.24) 
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and the leading asymptotic term Wo was computed explicitly in terms of the function 



It reads 



-^log(l+7i^(A)), i^± = iy{±q). 



(3.25) 



Wo(m, M) = - /"[imp^(A) + 5'(A)]i/(A)dA-J]{z.2log[msinh(2g)p^(a(?)] -logG(l,i/,)} 



1 fu '{x)u{^^)-l.{xy{^,) 

tanh(A — n) 



where G{1, z) = G{1 + z)G{l — z) , and G{z) is the Barnes function defined as the unique solution 
of the equation 



G{l + z) = Tiz)Giz), with G(l) = l and -^logG(z)>0, z>Q. 

az-^ 



(3.27) 



It was proved in that this leading term Wo gives rise to the leading term I^"* [Qs] of (I3.16P 
through the procedure described above. 

In [46] were also computed the first oscillating and non-oscillating corrections to (13.24(1 . Non- 
oscillating corrections to (|3.24l) lead to non-oscillating corrections in (|3.16ll . and we a priori do 
not need them for our purpose. We have seen however that, even if the first oscillating correction 
(|3.17p in (|3.16l) gives rise, through the process of summation that we will see in the next section, 
to a subleading term in the expansion for (e^2,nN,^ it 

may become leading once the difference- 
differential operator (|1.4I) is applied. Therefore we also recall the explicit expression of the first 
oscillating correction which was computed in [1^. It reads. 



Wi(m, [v]) e«™[po(g)-Po{-g)] + [j,]) g-im[po(g)-po(-g)] ^ 



with 



and 



W±i(m,M) 



= n 

(T=± 



v^v^ u{q) 



±1 



sinh {2q) p'^iq) p'q( 



m 



(3.28) 



(3.29) 



00-9(0-1?) 



[sinh(2g)po((Tg)] 



2u^ 



r(^^.) 



exp 



-2a 



KA) 



tanh((7 — aX) 



dX 



(3.30) 



The asymptotic estimates given above are uniform over 7 for 7 small enough: for all e > 0, there 
exists 7o such that, V7 < 70, 2|z/+ -|- < e and the corrections to (|3.28|) are of order less than 
0(m-3+=). Therefore (lOHIl indeed gives rise, through the process described above (see (I3.23P ). 
to the first oscillating correction Og (|3.17p of (|3.16|) . 
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Remark 3.2. We would like to draw the reader's attention to a remarkable relation between the 
terms Wo and W±i: 

W±i(m, [jy]) e^*™[^o*^^)~P'''^~^)] = e''^o'^™'[^^-'-l)^^o(m,[i/]) _ (3.31) 

Hence, the terms W±i partly restore the original periodicity ^ + n, n G Z, of det[/ + 7F], 
that is broken if one considers the leading term Wq alone. 

3.2.2 Leading asymptotic behavior of cycle integrals 

We have seen how to generate the leading term /i^'' [Qg] and its first oscillating correction Os [Gs] , 
in the asymptotic expansion (|3.16|1 of the cycle integral (|3.1I) : /i'^^[^s] is generated through 
(|3.23p by the leading asymptotic part Wq ()3.26|1 . whereas Os[^s] comes from (|3.28|1 . It remains, 
in order to be able to sum up their respective contributions in (|3.15l) . to describe more precisely 
their action functional of Gf 

For the purpose of the summation of the series (|3.15p . it will be convenient to distinguish 
the ^-dependent part of Wo, presenting the latter as a sum of two terms Wo = wj"^ + wj^\ 
where 

wS^\m, H) = Wo(m, M) and wJ^^M = - / g'iX) z^(A) dA. (3.32) 

9=0 J 

It will split accordingly the functional li^^ [Gs] into two parts: 

I^'>^[Gs] = {Hs + Ds)[Gs], (3.33) 

where -ffs[^s] is the part of /i''''[^s] issued from Wq^\ whereas -Ds[(/s] comes from W^^\ 

To obtain the explicit expression of Dg as a functional of Gs, one should proceed as follows: 
compute the s^^ 7-derivative of Wq^''[z/] at 7 = 0, 



q 



(9)r 



^g'{X)F%\), (3.34) 
7=0 J 2m 



(s-1)! ^ 

-q 

substitute in this expression F{X) and g{X) by (|3.22|) . 

i^5;w(^)M|^^^ = / ^.UA)[^.(A)e-W]% (3.35) 

-q 

and take finally the sum over k as in (|3.23|) . This gives 

(3.36) 



^sm = ]^^3.Gs^ 



{A}^ 
A + e,{Ar-i 



e=0 
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The notation {A}* means here that all the original cycle variables Ai, . . . , As are now set equal to 
the same value A. The variables Z2, . . . , -Zs are likewise equal to A. Notwithstanding, one should 
set here zi = A only after taking the derivative with respect to this variable. 

One can derive similarly explicit formulas for the action of Hg and Og. However, for general 
functions Qs, the corresponding expressions are rather cumbersome (see [46]). Therefore we 
restrict our analysis to a specific class of functions. Let us suppose from now on that Qs can be 
presented as 

.(i:j)^n^.(.,i:i)-..(i:j), 

where the functions 

cKiLl I^M, $2 I 1, (3.3^ 





are symmetric functions of the s variables {A} and of the s variables {z} separately, and satisfy 
the reduction property 



$1 u 




\ {z} \ Zj 




Note that, if the set {A} coincides with the set {z}, then the function ^2 becomes a constant 
whereas becomes a one-variable function: 



^2 ^ = ^2 = const, $1 cj I = [ = $1(0;). (3.40) 





Observe that the functions V and W in (|2.35p are precisely of this type and, hence, ^\j^^y\ is 
a particular case of (|3.37|) . Moreover, in the process of consecutive summation of the series 
(|3.15p . we will permanently deal with functions of the form (|3.37|1 . A remarkable property 
of such functions is that, regardless of their multi-variable nature, everything happens in the 
asymptotic regime as if they were a pure product of one variable functions. 

Let us give now the explicit form of the action of the functional Hg and Os on such functions. 



Proposition 3.2. Let Qs he a function of the form I13.37\) . Then 

■ "^2 , (3.41) 



7=0 



where 



z)(c^) = --^log[l + 7«>i(u;)], (3.42) 
and with ^i{uj) and ^2 defined as in (iS.40\ )- 
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Proof — The explicit form of the functional Wq'^^[i^] follows from (|3.26l) . However, to prove 
Proposition 13.21 it is enough that it can be written in the following quite general form: 



n=o j=i 



(3.43) 



where z^(^) is given by (|3.25p and T„ are some functions or distributions. In order to obtain the 
corresponding part of the asymptotics of the cycle integral of length s, one should first take the 
s^^ -(/-derivative 



E E' n%r|,.„/r»(«..-.«")n^''fe)<i"« 

n=0 ki,...,k„=lj=l ^ _ j=l 



7=0 ^— ' ^ ^ kA 

n=0 ki,...,k„=lj=l 



(3.44) 



where vq is given by 



-1 



Ztti 



(3.45) 



and prime means that the summation over /ci, . . . , /c„ is taken under the constraint 'Y^=i ~ ^■ 
Then, using (lOTTl and dH^S]), we obtain 



(-1) 



o_i CO s n nki 



is 



— V V' TT ' 
1)1 ii t.i 

' n=0 ki,...,kn=lj=l ■' 



9 ^ 



It remains to observe that, due to the reduction property (|3.39|1 . 



\ {6}'^S...,{en}'="^ 

Comparing with (|3.44|) we see that it gives (I3.41|) with z> defined in (I3.42p . 



{in} 



k id"e- 



(3.46) 



□ 



Proposition 3.3. Let Qs he as in (fO^). Then 
(s - 1)! 



7=0 \ ±g 



(3.47) 



where 



pW(a;) = -— lo, 
27r« 



1 



1 + 7^>i u 



(3.48) 
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Proof — The proof is very similar to the one of Proposition 13.21 Consider, for example, the 
oscillating functional O^. The action of this functional is generated by the oscillating term 
W+iM e*'"[Po(9)-Po(-'?)]^ which in its turn is expressed in terms of W^^V] by dOTjl : 



Once again, we do not need the complete explicit expression of Wq^^ ■ Let simply 
r„_ii_i/i;(0)r 



(3.49) 



(3.50) 



in which we have stressed that the oscillating term of the asymptotics of the Predholm deter- 
minant is proportional to ^{q) ^{—q). For the remaining part T[h>], we can use a representation 
similar to (13.4311: 



«=0 j=l 

where T„ are some functions or distributions. The s^^ 7-derivative then reads 
d;{uiq)ui-q)T[u] ' " - 



(3.51) 



n=0 k-,ki=l 



7=0 



/(I. 
r„(6,...,en)F'+('z)i^'-(-g) n^''(^^) (3.52) 

where z^o is defined in (|3.45|1 . and prime means that the sum over ki,. . . ,kn is taken 

under the constraint + A;_ + X]j=i = •^^ Then, using (|3.2ip . (|3.22|) . and multiplying (I3.52p 

by e9(i)-9(-q)^ we obtain 



00 s rJt+ rsk- n „kj 



k+\ k.\ 

X j T„(ei,...,en) Qs 



k,\ 



7=0 



d"^ (3.53) 



It remains to observe that, due to the reduction property (|3.39|) . 

\ q, iqr^,i-qr- m 
I -q 



-q, M'=^{-<?f-^{el}'^•••,{en}'=" 
q, {(?f^{-<?}'=-^{el}'^•••,{en}'=" 



(3.54) 



Comparing with (|3.52|) we arrive at the action (|3.47l) with z>(^) defined in (I3.48p . 



□ 
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Thus, the action of the functionals Ds, Hg and Og is exphcitly defined at least on the class 
of functions Gs admitting the representation (|3.37p . In that way we have found the leading 
non-oscillating and oscillating asymptotics of cycle integrals. However, this is not yet enough 
to obtain a correct estimate of the infinite series of multi-cycle integrals (|3.15p . The matter 
is that some a priori sub-leading corrections after their summation may give non-vanishing 
leading contribution. We will consider this question in more details in Section 14.21 where the 
mechanism of summation will become clear. Here we merely would like to note that some 
part of the corrections could easily be included into the action of the functionals considered 
above. For example, we did not use the explicit form of the functional Wq^^ [u] in the proof of 
Proposition 13.21 Hence, instead of >Vq°^[z^] alone, we could as well have considered all the terms 
in the asymptotic expansion of log det [/ -|- 7F] which can be written in the form (|3.43|1 . and the 
statement of this proposition would still be valid. 



4 Asymptotic summation of cycle integrals 

In this section we sum up the series (|3.15p in the asymptotic regime m — 00. We already labeled 
each pure cycle integral I(s,p) with respect to the position of the variables on which it acts. We 
define similarly -D(s p) (resp. -ff(s,p) and as the operator that acts on the variables of the 

p^^ cycle of length s as Ds (resp. Hg and Og), and by doing nothing to all the other variables. 
Then, if we also define Rg to be the functional given by the remaining part of (I3.16P (which is 
in principle of order o(l) for non-oscillating corrections, and of order o(^) for oscillating ones), 
and R(s,p) to be its conterpart acting on the cycle (s,p), we can reexpress (|3.15|) as 

Gi...n(7) = E n -rr n i^is,p) + H^s,p) + + Ris,p)] * ^\J^,y\ , (4.1) 

£l,.../„=0 s=l p=l 
/ 1 ^s — i^s 

where Ug = - — — — —. There, just as in the case of multi-cycle integrals, the order of the 
different operators is not important due to the symmetry of the function in its two sets 

of arguments {A} and {z}. 

Before starting to sum up explicitely the series (14.11) . let us make a usefull remark: since we 
only need the n**^ 7-derivative of Gi,,,n at 7 = (see (|3.14p ). it is enough for our purpose to 
obtain a ^-equivalent form of the result. 

Definition 4.1. Two functions v?i(7) and v?2(7) are said to he 7-equivalent of order n if 

fc = 0, l,...,n. (4.2) 



a>i(7) =5^^2(7) 
7=0 



7=0 



Therefore, at any stage of the computation, we may replace Gi,„n by some 7-equivalent 
function of order n which, by abuse of notations, will still be called 
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4.1 Summation of the action of H, 



Let us start by summing up the action of the operators -f^{s,p) • We first expand the products of 
the operators in (|4.ip according to the binomial formula: 



1 



E 



1 



rj(4-r.)! 



n [-^(s:p) + ^(s:p) + ^(s,p)\ n -^(s,p) • (4-3) 



As we already mentionned, such an operation is possible due to the symmetry property of the 
function -^j^^ji on which the operators in (|4.3p act. Substituting this into (14. ip and changing 
the order of summation, we arrive at 



Gi...„(7) 



E n 

ri,..., r„=0 s=l 

£l,...,£„=OS = l S = lp: 



n n [^(^-p) + ^(^.p) + ^(^-p)] 

s=l p=l 



Here, for convenience, we have divided the set of arguments of into two subsets: the 

operators and R(s,p) act on the variables {Ajj^^j and {z}j^^y, while the operators 

i?(5,p) act on {/i} and {y} j^^j . 

Recall that each individual -ff(s,p) acts non trivially (as the functional Hg) only on variables 
IJ's,p,j and ys,p,j, with j = 1, . . . ,s. Therefore one can write 



{fJ-s,p,j}i<j<s U {A}j, , U {/i},srp 

kMl+l^wl 1 {y,,p,,}i<,.<, U {z}j,. U {yk.rp 



(4.5) 



where Hg acts on the first sets of variables {/Us,p,j}i<j<s and {ys,p,j}i<j<s according to (|3.41ll . 
while all other variables can be considered as auxiliary parameters. The set J^^ is defined as 



(4.6) 



It is clear that j^^jl+j j^^^l has the form (|3.37ll . with 



Vi 



J{r}\ + \J{1}\ 



UJ 



{/^S,p,j}l<J<f 

{y s ,p ,j} i<j <s 



(4.7) 



$9 



n 



Vi 



k{r}l + l>^{£}l 



UJ 



{f^s,p,j} 



i<i<s 



{/^S,p,jr}l 



i^eiAjj U{m} rj, 
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where, for simplicity, we have omitted the additional sets {A} j, , U {^} ,srp and {z} j, , U {y} ,175, 

-'{1} ' -^{1} 

in all the arguments in (|4.7|) . (|4.8p . Applying Proposition 13.21 we obtain 

(-l)-i 



H{s,p)*^\j{,}\+\j{t}\ 



where 



z.«'PH = — log 

ZTTl 



I + 7V, 



7=0 [-'{rjl+IJ^;,!! 



"'{r}l + l^{f}l 



{1} 



(4.9) 



(4.10) 



The obtained result still admitting the representation (|3.37p . we can act with all operators 
H(^s,p) consecutively. It gives 



n is 



nn^(^.p) *-^-^{.-}i+i^wi 



n 



s=l p=l 



7=0 



^1 



where 



z)(^) = ^lo, 



Substituting this into (|4.4|) and summing up over all is we arrive at 

00 n / s„ n rs 



Gi...n(7) = E n n n + ^(..p) + R(s,p)] * 



ri,..., r„=Q s=l 



where 



(H) 



\-J{r}\ 



S = \ p=l 



. S = l 



(0), 



7=0 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



Since it will turn out to be rather important for the summation of the remaing terms, we have 
here explicitly indicated that the functional yV^^^ {m^ [i>]) depends on the distance m. 

We can now perform the summation over s in (|4.14p by extending it up to infinity, which 
means that we replace T^^^ by the following 7-equivalent function, still denoted by T^^^ : 



\J{r}\ 



(4.15) 



Thus, the successive action of operators -f^(s,p) produces a complete decoupling of parts of 
the variables, and eventually restores some of the Fredholm determinant asymptotics (|3.26l) . 

Remark 4.1. Note that the summation performed here is exact, and that we did not use the 
explicit expression of >VQ'^^(m, [i)]). We could as well have considered the action of some operator 
Hg including not only the leading contribution Hg, but also subleading (non-oscillating) ones, 
provided they originate from the same kind of terms (|3.43ll . In fact, the very same process will 
enable us, in the next subsection, to sum up some part of Rg as well. 
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4.2 Contributions from R(s,p) 



Unlike the original function (I2.35p . the obtained function -^j^ j| depends on the distance 

m. Let us specify this dependence. We have 

q 

Wj°)(m, [u]) = -im I p'o(A) i>(A) dA - ^ z)^ log [m smh{2q) p'^iaq)] + C[u] , (4.16) 



with 



= 2 J tlliX- J) ''^^ + W tanh(., - A) ^^ + 
in which vioj) is given by KW\ . Then the function J^y^ | reads 

J^lf^l^^ = exp I -im J p'oiio) uiio) dtjl • JJ [m sinh(2g) Poiaq)]"^^^""'^ ■ e^'^l 



(4.17) 



-q 



<T=± 



In order to complete the calculations, we should in principle sum up the action of the operators 

( H) 

£)(s^p) and 0(s,p) on ^| j | , provided that R(s,p) produces only suhleading corrections. This was 

the case when acting on m-independent functions of the form (|3.37|) but, due to the fact that 

m, tl 
.{H) 

\J{r)\ 



( H) 

T\j I now depends on m, this is no longer true. Let us explain why. 

The new function lF}j ^ , is still of the form (I3.37p . but now $2 depends on m as 





= exp |m $ ^ j + logm $ ^ I }> • «>2 I )l\ ) , (4.19) 

with $, $ and $2 satisfying the reduction property (I3.39p . We see that the differential operator 
in -0(5 p), when acting on such m-dependent functions, may produce terms that are linear over 
the distance m: hence, although this operator originates from a term of order 0(1) in (|3.16l) . it 
behaves as a term of order m when acting on (I4.18P . The same phenomenon may happen for 
the action of R(s,p)- if the latter contains some derivative operators, then these operators may 
produce non vanishing contributions to the final answer although they are effectively subleading 
when acting on each cycle; for instance, a differential operator gives a term of order O(^) 
when acting on the original function |, but it gives a contribution of order 0(1) when 

acting on the new function ^ ^ | which is the result of summation over all cycles of the action 
of on Therefore, we need more informations on the structure of the subleading 
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terms in (|3.16|) in order to be able to extract from Rg the part of the corrections that eventually 
contributes to the leading order. 

The structure of the series (|3.16|) for was studied from the asymptotic expansion of the 

Predholm determinant (|3.19ll and related Riemann-Hilbert problem in [46]. This analysis being 
extremely cumbersome, we do not have the explicit expression of the series (j3.16ll . However, 
we obtained in |46j some information on the generic form of the corrections: they are given 
in terms of the function Qs or of its derivatives, the latter being evaluated at the endpoints 
±q or integrated, together with some weight function, on the interval [—q,q]. Therefore the 
functional Rg indeed contains differential operators. However, it was proved in [1^ that the 
non-oscillating correction li"' "^"^^^ [Qs] of order n in ()3.16|1 contains at most derivatives of order 
n of Qs^ whereas the oscillating correction li'^'°^^\Qs] of order n contains at most derivatives 
of order n — 2. They correspond respectively to the contribution of some differential operator 
_2_^{",nosc) jjj ji^^ ^(n,riosc) j^gjj^g q£ Qj-^gj. most n and 5}^"'°'^'^^ of order at 

most n — 2. Therefore, only the part of maximal order of these operators may produce, when 
acting on functions of the form (I3.37|) . (I4.19P (and more precisely on the exponent), contributions 
on the same order as Dg or Og. These are the contributions that we need to determine. 

For convenience, in order to distinguish the effect of non-oscillating and oscillating correc- 
tions, we define two functionals Rg°^'^ and R°^^: R'^"^^ contains all the non-oscillating corrections 
with respect to and R°^^ contains all the oscillating corrections with respect to Og, in the 
sense of the asymptotic series (l3.16|l FI. 



4.2.1 Non-oscillating contribution 

Let us first discuss the effect of the non-oscillating part Rg°^^ of the corrections. Clearly, as we 
already noticed, at each order n, only the differential operators ^2)1"''^°'^^'', and more precisely 
the part of which results in the maximal order n of derivatives, may 

eventually produce, when acting on m- dependent functions (|4.19p . contributions of order 0(1). 
For n > 1, the explicit expression of the differential operator Xi^'^"^^^ of maximal order n is not 
known, but it can be shown from [l6] that it has the form 

s 

J2ikj+(j)=n, (4.20) 

^i=^3='"l j=l 

where Ca{{k}, {£}) are some computable coefficients. Here the action of the derivatives (on 
some function Gg) is followed by the evaluation of all the variables z and A at the same point 
aq, a = ±. 

Obviously, when acting on some m-dependent function Gg of the type (|3.37|) . (I4.19p . this 
operator can produce 0(1) contribution only if it acts on the part of the exponent which is 

^It means that we choose some arbitrary integer no and consider the asymptotic development for and 

RT'' up to O(log"m/m"o+i). 
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linear in m, i.e. on 
exp < m <I> 




(4.21) 



Using the reduction property (I3.39|) for $i and $2 and the fact that the contributions of order 
0(1) can be obtained only in case each derivative hits exactly once the exponent (|4.2ip . we 
obtain 



a=± 



Y.'CA{k},m)ll[9'i^^P[-9'i>^j)]' 

{k},{e} j=i 



1 + 



log m 



m 



(4.22) 



n. Here we 



where the prime means that the sum is taken under the constraint Yli^j + 
have used the conventions (|3.40|1 for the reduced functions '^i and $2, and we have defined the 
one- variable function g' to be given in terms of the reduced function <I> as 



g'iX) = 9, CD 



A 

A + e 



(4.23) 



e=0 



Hence, everything happens (at least as far as the leading term is concerned) as if we would be 
acting on functions Gs of the type 




$2 • n ^i(^i) • exp < m ^ [gizj) - giXj)] 



(4.24) 



Thus, the problem of finding the structure of the 0(1) contributions coming from the action 
of Rg°^^ on some m-dependent function Gg of the form (|4.19p reduces to the computation of the 
leading behavior of the functional action of Is on Gs (I4.24P . Strictly speaking, the functional 
argument of the cycle integral 2^s[^s] should not depend on the distance m but, in the case 
(I4.24p . we can simply re-define Po ^ p = Po — ig and Gs ^ Gs = Y\j=i (-^i) (^^ recall that 
$2 is just a constant that will be in factor of the result). This is then a special case of (|3.18ll . 
and therefore the leading contribution of the corresponding cycle integral can be obtained from 
the s^^ 7-derivative at 7 = of 



-im 



po- 



/ p'(A)z>(A)dA-^{i>2log[msinh(2g)p'(ag)]} + C[z>], (4.25) 



where u is given in terms of $1 as in (|3.42l) . The fisrt term of the r.h.s. of (I4.25p . which is 
linear in m, corresponds to the action on (|4.24p of the operator Dg that we already know. It 
is the other (constant in m) p'-dependent part that eventually reproduces the contribution we 
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want to compute, namely the 0(1) contribution coming from the action of the whole series 



X^n>i ^2)^"'''^°'^'^^ on m-dependent functions of the form (|4.24p . It reads 



^2, 



(4.26) 



7=0 



where we recall that £> is given in terms of as in (|3.42|) . and that 

A 



p'{X)=p'o{X)-id,<!> 



A + e 



(4.27) 



e=0 



We have thus determined the explicit contribution of order 0(1) coming from the action 
of Rg"^'^ on m-dependent function Qs of the type (|3.37|) . (|4.19|) . Let us call this quantity 



^(1, nose) g3jp}j(,j|;g}y^ Rs'''^"^^'[Gs] is given by the expression (|4.26p . If moreover we denote 



,(1; nose) r 



byi? 



(sub; nose) 



Gs] the remaining part of Rg°^^[Qs], we have 



R^°''[g,] = Ri^'^'^'^^igs] + R'i'''^'''°''^[Gs]- 



(4.28) 



Remark 4.2. What we have done here is to act formally with the asymptotic series (|3.16l) on 
m-dependent functions Qs of the form (I3.37p . (I4.19p . and to identify, at each formal order n 



in m, the part of zzwl: 



1 j{n; nose) 



[Qs] contributing to the order 0(1). We know that the summed 



contribution of these leading parts, that we have called R, 

1 j{n; nose) r 



(1; nose) 



Qs], is given by We 

also know that the remaining terms in each -^Is"''^"^^\Qs] are of order 0( '°f7" ). However, it 



is not obvious to prove rigorously that the sum 

Rr''[Qs] of these remaining terms stay 
subdominant in the final answer. It should follow from the analysis presented in [3^ but here 
we just assume that it is so. 



It is now possible to sum up the successive actions of the corresponding operators R 



on the function T^"'^ KTEtl . The action of R', 



(1; nose) 



(1; nose) 



being similar to the one of Hs (compare 



(I4.26P to (|3.41|) ). we can apply exactly the same procedure as in the previous subsection, and 
the action of R^(l'p°^'^^ exponentiates similarly as in (|4.14|) . Therefore, a 7-equivalent form of the 
series (|4.13p is 



Gi...„(7)= 74n[^(«.p)+0(^.p)+ 

where 



ri,...,r„=0 p=l 



j~,(sub; nose) pose 



,( ), (4.29) 



jr{H) ( {A}j{,j 



n 



(4.30) 
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Using the explicit expression (lil8D for we obtain for the new function J^^^^: 



(R) 



\J{r}\ 



exp ^ —im J p'q[u)) u[u) 



duj > • [m sinh(2g) p' (erg)] 



X Wi 



J{r}\ 



{r} 



n V, 



{r} 



in which z> is still given by (I4.12|) and 

g 



H + e, {z}j^^^ 



duj . 



e=0 



(4.31) 



(4.32) 



4.2.2 Oscillating contribution 

Let us also compute the leading action of the oscillating part R°^'^ on m-dependent function Qs 
of the type (|3.37p . (I4.19p . Similar considerations based on the explicit structure of the oscillating 
corrections in the asymptotic series (I3.16P [i6] allow us to decompose Rf^[Gs] as 



(4.33) 



Here [Qs] takes into account the contributions to the leading oscillating order of the action 

of each subleading oscillating functional in (|3.16|) . and therefore is an oscillating contribution 
of the same order as Os[^s]) whereas the remainder i^i^"*^' [0^] is subleading with respect to 
Os[Qs]- Like in the previous case, the modification at the leading order reduces to the extension 
of the term logpg to logp' = log [p'q — j^d^ . Namely, let us define Og = Os+R^}'^^'^'' = 0^+0~ . 
Then, for Qs being given by (I3.37p . (|4.19p . we have 



JJ [msinh(2g)p'±(c7g)] 



±2!>(±)((Tg)-l 



where is still given by (|3.48p and 



±g + e± 



(4.35) 



e=0 



The successive action of the oscillating operators 0(s,p) will be summed up perturbatively in 
Section KM 
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4.3 The action of -D(s,p) as a continuous generalization of multiple Lagrange 
series 

Let us now sum up the successive action of operators -D(s,p)- Applying again the binomial 
formula, we obtain 



oo n 



ri,..., rn=0 s=l s=lp=l 



X E n ^ n ft °<... • ^'i:;.H.„„ ( j*;'; ^ ) • 



in which the operators -D(gp) act on the variables {A} and {z}, whereas 0{s,p) and R[s,p) 
R^^^^^^°^^^ + R^^^p^^^^"* act on {ji} and {y}. We therefore need to compute 



^(D) / {lAj(r} \ V- TT (^s) 



,...,£„=Os=l 



This series is actually a continuous generalization of the multiple Lagrange series. 
Recall that the standard Lagrange series has the form (see e.g. ^81j ) 



°° 1 d" 

n=0 



(4.38) 



e=0 



where F{z) and 0(2) are some functions holomorphic in a vicinity of the origin. If the series 
(|4.38p is convergent, then it can be summed up in terms of the solution of the equation 



z - (t){z) = 0, (4.39) 
and the sum is given by 

It is possible to generalize the Lagrange series (14.381) to the case of functions (j) and F 
depending on several variables, and even to consider the corresponding continuous limit. The 
sum of such generalized series is then expressed in terms of a solution of an integral equation. 
All the details about these generalizations are given in Appendix O We just apply here the 
result of this appendix to our particular case. 
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Substituting the action (I3.36P into (|4.37|) and setting k = \J{r}\, we obtain 



^ 1 



(D) ( {lA 

\ {y} 



oo n 

E n 

£i,...,£„=0 s=l 



2m 



XT 



(R) 



I {y}, U ( {Xs,p + e,,p} U {A,,p} 



s-1 



• (4.41) 



Recall that the notation {As^p}* means that the variable A<j,p is repeated s times. Moreover, here 
and in the following, we use simplified notations for 

n 

U- U . E-EE- (412) 



s,p l<s<n 
l<p<ts 



s;p s=l p=l 



Using now the explicit expression ()4.31|) of \ i we have 



E ns/{nn 



oo n 



% ( n Is 



ei,...,en=o s=i 



dA 



=ip=i 



■ d 

27ri 









1 


exp < 


— im y 




es,p=0 





JJ [msinh(2g)p'((j(7)] 



Wi 



e\+k 



<T=± 

k 

t=l 



nn^i 



u {A„,4 

u {A„,4 

A,, 



U,V I S— ip— i \ 



, (4.43) 



where \i\ = YTs=\^!i- -^^ order to lighten notations, we did not write explicitly the arguments 
{/i} and {y} of the functions V and W . We also have omitted most arguments of the functions 
V and -p' . For instance. 



V\ijJ ] = V \ ijj 



27ri 



log 



1 + 7V w 



{/i}, U {\u^v} 

u,v 

{y}, U {Xu,v + eu,v} 



(4.44) 



and one should understand p'{crq) in the similar way. 

Observe that we have, at most, to compute the first order derivative with respect to each 
variable e^^p. It is thus enough to linearize the arguments of the functions V and W in the 
vicinities of e^^p = 0. The linearized version of the function V reads 



Vi 



e\+k ^ 



U {A.,p} 
U {As p + ^s,p} 

s,p 



exp 



h-iJ2 ^«'P ^ - ^-'P) + ^k{uj)\-l, (4.45) 

I s,p ) 
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where K{uj — A) is given by p.lOj) . and the function ^'^(w) by 

{/^} 



{y] 



s\nh.{uj — fit + iC) sinh(a; — yt — iC) 



k ^' 



sinh((j -yt + iC) sinh(w -fit- <) ^^-^ 



The hnearized form of 0{uj) follows immediately from (I4.45P : 



^ / K{uj-X)dX. 



(4.46) 



yt 



exp /? - i ^ es,p K{uj - A^.p) + ^fcl^^) - 1 



(4.47) 



As for the function W , its linearized form is rather cumbersome. On the other hand this function 
will not play a significant role in this section, therefore we give here the linearization of in a 
rather formal way, namely 



U {A.,p} 

s,p 



s,p 



^t,r 9 (As,p, At^r); • • • 1 ) (4.48) 

\ s,p s,p t,r / 



where g^^\ g^'^\ . . . are some functions of Xg^p. The explicit formulas for them will be given in 
Section 15.11 Then we have 



oo n 



% ( n Is 



E n^/ nn 



fi,...,^„=Os=l 



: JJ [m sinh(2g)j5' (erg)] 
k / 

n/i 

t=i 



s=l p=l 



exp < —im 



1 



u> ) duj 



n ts 



/^t I X] '^u,v K{flt - Xu,v) 1 W W ^Z'* ( ^s,p I -i ^ iu,v K{Xs,p - Xu,v) I , 
V u,v / s = lp=l \ u,v / 



where the notation {{g^^^}) is a compact form of ()4.48ll and where we have set 
/(A|2;) = e^+^+*^W-l. 



(4.49) 



(4.50) 



The representation (|4.49|) is exactly the continuous generalization of the multiple Lagrange 
series considered in Appendix O (see (IC.19P ). Thus, we can simply apply the results derived in 
this appendix to the concrete case of functions /, Ffc- 



43 



Let us give the final result: a 7-equivalent form of the new function J'f'' dOOll reads 



im 



p'q{ijj) z{uj) duj > [m sinh(2g)p' (erg) 



-z2(o-g) ^ c[-z\ 



t=i V -<? / 



det 



/ + iK{i - A) U I / K{u; - z{io) du; 



Here the function z{X) solves a non-linear integral equation, 



(4.51) 



^(A) = hlX\-i K{LO-X) z{uj) duo , 



(4.52) 



with /s given by 



/s(A|x) = — log 



1+76 



(4.53) 



The symbol means the derivative of /s(A|x) over x. The Fredholm determinant in the 
denominator of (I4.5ip is the Jacobian of the equation (|4.52p (compare with (|4.39|1 and (|4.40p ). 
The functior0p'(A), which appears in (|4.51|) taken in the points A = ±q, satisfies a linear integral 
equation, 



p'{X)+ / Kj:iX,co)p'iu;)du;=p'o{X) 



(4.54) 
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where 



/..(A,.)^.A-(A-.)/iL|-./A-K-.).,Od« 



(4.55) 



We would like finally to mention that the functions z{X) and p'{X) depend on the parameters 
{//} and {y} through (see (|4.53ll ). 



z{X) = z{X 



{y} r 



p'{X)=p'[X 



M r 



(4.56) 



and still satisfy the reduction property (|3.39l) 



^ We use the notation p'(A) for this function, since eventually p(A) will have the sense of the dressed momentum 
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Let us now explain how to obtain this result from the results of Appendix O 
According to Appendix IC.41 the final answer for is given in terms of a function z^"^ (A) 

solving an integral equation (see (IC.25P ) 

z(")(A)-/E„|^A|-z|i^(..-A)zW(a;)d^j =0, /^JA ] x) = |^ ^ (A | x). (4.57) 

Following our usual strategy, we can replace this equation with a 7-equivalent one by sending 
n ^ 00 in /s„ as Us gc 7*. Then the function z^'^\n) is replaced by z{^) satisfying the integral 
equation (|4.52|) . since 



00 / i\S-l 



lim /s„(A|x) = V 

n— >oo '—^ 

s=l 

1 



,/3+a:+*fe(A) 



1 



27118 

log [l + 7/(A|x)] =/s(A|x). 



(4.58) 



Comparing (I4.5ip and (I4.49p one can observe that the action of the operators ^(s,p) leads to 
the replacement of the functional argument z>(A) by —z{X). This happens due to the fact that 
the linearized form of v (|4.47p coincides with the r.h.s. of the equation (I4.52|) up to the sign 



U{A.,p} 



u to 



s,p 



U {Xs,p + es, 

s,p 



(4.59) 



Since the summation of the continuous generalization of the Lagrange series reduces to the 
replacement of the sum over es,p by the integral with z{X) (see (|C.24j) ). we obtain 



z>(w) ^ -/e cj I -i / K{oj - A) z{\) dA = -z{uj). 



(4.60) 
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Finally, let us explain how the function p'[X) and the equation (|4.54p arise. From (|4.32|) and 
(Onll . we havcil 



1 

p'(A) ^p'(A) =p'o(A) + j M d,z(^u; 



A 

A + e 



(4.61) 



e=0 



Differentiating the integral equation (I4.52p we find 



dfZ \ LO 



A 

A + e 



A 

A + e 



de = -i^s(A,c^) 



(4.62) 



e=0 



'Here also, we do not specify that z depends on the extra sets of parameters {/^} and {y}. 
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Here —K-s{X,lo) in the r.h.s. appears due to the derivative of the function ^'^ (I4.46p . Hence the 
derivative d^z can be expressed in terms of the resolvent R^{X,u)) of the integral operator K^: 



dfZ \ LO 



A + e 



e=0 



Substituting this into (I4.6ip we arrive at 



(4.63) 



(4.64) 



which means that p'{X) solves the equation (|4.54p . 



4.4 Oscillating corrections 

The series (I4.36P takes now the form 

oo n ^ \y. n 

oi...«(7)= E n^nnt 



o 



+ R 



=0 s=l 



with ^1 ^ I given by (I43TD . Note that J^^j' , is still an m-dependent function of the type (|3.37|) . 



s=l p=l 



-(D) 



(D) 

J{r}\ 




(4.65) 



l'^{r}l 



(|4.19p . and therefore we can use the results of Section 14.2.21 

At this stage, the functionals which remain to be summed up produce only subleading 
contributions. Hence, if we were only interested in the leading term of (e^^™), it would be 
enough to consider only the leading term of these series, namely the term ri = • • • = r„ = 0. 
Recall however that we need to keep the leading oscillating contribution since, after second 
lattice derivative, it may produce a term competing with the leading non-oscillating one. We can 
nevertheless forget completely the remaining terms R(s,p) > since they merely produce corrections 
to one of these two main contributions. Moreover, since the oscillating operators 0(s,p) are 
themselves of sub-leading order, it is enough to restrict our analysis to the terms of (I4.65P that 
are linear over 0(s,p), namely the terms for which one takes the value or 1 while all other 
are equal to 0. We thus obtain 



Gi...„(7)=-^; 



(D) 



+ y^usOs 




+ corrections, 



(4.66) 



in which the corrections are either non-oscillating corrections to the first term, or oscillating 
corrections to the second term (note that the non-oscillating corrections may be of order higher 
than this second term, but they are nevertheless not important for the final result). Observe 
once again that, due to the fact that cx 7'^, the summation over s in ()4.66l) can be extended 
to infinity. This will be done in the end of our calculations. 
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To complete the summation of the series, we should now calculate the functional Og [J^s^^] 
Let us present Gi...„(7) in (I4.66P as the sum of three parts 



corrections . 



(4.67) 



a=0,±l 



where 



fc = l 



■J 1 



(D) ( {^} 

I {y} 



(4.68) 



The non-oscillating part G^^\{'y) has in fact already been computed. It is given by (|4.5ip 
in which one should set {^} = {y} = 0, which is easy since we almost did not specify this 
dependence: basically, the parameters {/x} and {y} enter the functions z{\) and p'{X) (I4.56P 
through the function ^'^ (I4.46j) . Thus, in order to obtain Gf^\{'y), it is enough to set ^'^ = in 
(ICTIl and A; = in (li^TTl . 

Let us now compute the oscillating corrections g[^^„(7). The new function J'jf^ KETh is of 
the type (jOTll . KWh . with 



$1 e 



{^} 

{y} 



(4.69) 



where /(A|a;) is given by (|4.50|) and 



{y} 



i / Pq{u>) z{u)) duj 



(4.70) 



The functionals (|4.34|) send /ii = ^q, yi = ±g and = yj for j > 1. This means that the 
function z{X) turns into ^'^^•'(A) satisfying the equation 



z(±)(A) = f^'^ { X \ -i I K{oj- X) zW(a;) duj | , 
with given by 



(4.71) 



#Va I x) = — log 



1 + 76' 



,/3+x+1'(±)(A) 



and 



^r(±)(cj) = log 



sinh(w — fj. + i() sinh(w — y — K) 



sinh(u; — y + iCj sinh(u; — fi — i() 



y=±q 



1 

J K{uj - X) 



dA. 



(4.72) 



(4.73) 
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Similarly the function p'{X) (|4.54|) turns into p±{X) satisfying 



-9 



where 



4±)(A,a;)=iK(A-u;) J K - u) z^^\0 



(4.74) 



(4.75) 



Comparing (147721) with K4S\\ for i>(±)(A) we conclude that z>(=^)(A) = -z(^)(A). Thus, we obtain 



[msinh(2g)p^ (erg)] 



7=0 \ ±q 



(4.76) 



Finally, similar arguments as in Section 143) (see (I4.62p to (I4.64p ) enable us to show that the 
quantity p'±{(Tq) defined by (I4.35p . (|4.70|) is equal to p'_^{aq), where p±{X) is the solution of the 
equation (14.741) . 

Observe that the k^^ 7-derivative acts on an expression which is independent on k. Therefore 
the sum over k in (|4.68|) . once extended to infinity, gives the Taylor series 

00 

u, 0± [^f = exp {±im{p',{q) - p',{-q)) + C[-z^^^ T 1] - C[-z(^^]} 



k=l 



X 



[m sinh(2(7) p'^((T(7)] 



a=± 

Combining now all the above results, we obtain 



gS11(7) = exp I im J p'^iX) {z(^\X) ± l) dA I J] ["ismh(2g)p'±(a<?)] 



-(z(±)((Tg)±l) 



/ si'^ (A) (A) dA; ... e^[-.(±)Ti] 



det 
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. (4.78) 



Here in complete analogy with (|4.48|) the functions F^^^ and {g±^} are defined from the linearized 
form of 

Tq , U {As,p} 



s,p 



±q , U {Xs,p + es^p} 

s,p 



\ s,p 



S,pj7 ■ ■ ■ 



(4.79) 
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5 Leading asymptotic behavior of correlation functions 



In the first part of this section, we gather all the previous results for the asymptotic behavior 
of (e'^^"). This will enable us, in a second part, to obtain the leading asymptotic behavior at 
large distance of the spin-spin correlation function (afa^^i). 



5.1 Asymptotic expansion for (e^^™) 



Let us indicate explicitly that the functions G^^\{'y) (|4.68p depend on the distance m and 
parameter /5: ^^'^^^(7) = G^^\{'y\(3,m). The remarkable property of the obtained 7-equivalent 
results is that they do not depend on n: GS^^„(7|/3,m) = G(<")(7|/3,m). Therefore the series 
(I3.14P is the Taylor series of g^'^HtIA 

m) at 7 = 1, hence leading to : 



= 5] G('^)(/3,m)[l + o(l)], with G^''HPM = fy_!^"!, - (5.1) 

The term G^°\p, m) in (15.1(1 gives the leading non-oscillating asymptotic behavior of (e^^™), 
while &^\p,m) describe its leading oscillating correction. Let us consider these two terms 
separately. 

5.1.1 Leading non-oscillating behavior of (e'^®™) 

It happens that, at 7 = 1, the non-linear integral equation for the function z{X) ^^5^ degener- 
ates into a linear one (taking also into account that, in this case, = 0): 

-q 

Setting z{\) = one gets the integral equation (il-14\ ) for the dressed charge: 

z{X) + J '^K{\-^,)z{^i) = l. (5.3) 

The linear equation (I4.54p also simplifies due to the fact that = 1/2-Ki at 7 = 1, 

p'(m) + ^ 1 (m - A) p'{\) dA = p'o(m) . (5.4) 

Comparing the last equation with (II. lip , we find that p'{\) = 2'Kp{X) and, hence, the function 
p{X) can be identified with the dressed momentum (see (11.12(1 ). 



-|- corrections 

7=0 
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Due to the same property of the Predholm determinant in the denominator of (I4.5ip 
turns into det[/+ ^K]. Thus, we have 



,,,, F„[£-,Mzip)^,:... 

[27Tp{q)m sinh(2g)] '^e'^ l^Tvii , ,2rr . ^ r.i ' (5-5) 



det^ [/ + ^K] 

where Z = Z{±q), and we have used that Z{X) and p{X) are even functions. 

It remains to pass from the symbohc form of Fq to a more specific one. When hnearized, 
the product in the representation (|2.20l) for the function W\£\ turns into 



nn 



sinh(2;s,p - A^/y - <) sinh(As'y - Zs,p - iC) 



i,p s',p' 



)\ sinh(zs_p - Zs',p' - iC) sinh(A5,p - A^/y - iC) 



exp < 



^s,p ^s',p' 5'"^'* {^s,p, As'y 

s,P s',p' 



where 



g^^\X,fi) = -- [sinh-2(A-;U + iC) + sinh-2(A-M-<)] . 
After the summation of the continuous Lagrange series, this part becomes 



nn 



sinh(z5^p — Xg'^p' — iC) sinh(A 



s',p' ^s,p 



^\ sinh(2^,p - Zs',p' - iC) sinh(As,p - A^'y - iC) 



■Co 



where 



Z(A)Z(/i) 



sinh {X — p — iC) 



dX d/i . 



(5.6) 



(5.7) 



(5.8) 



(5.9) 



The remaining part of W contains the Fredholm determinants (12.16(1 . whose kernels depend 
on the following products (see (PTZl) . (ITTH]) ) 



p„=n 



sinh(t(; — Xs^p + iaC) 



sinh(ti; — Zg^p + iaC) 



expiy^e,,p(7^^'")(A,,p)L a = 0,±l, 



s,p 



(5.10) 



where g^^'^'^X) = coth(i(; — A + ia(). Let us introduce the m-periodic Cauchy transform of the 
dressed charge Z, 

q 

(5.11) 

-9 



z(u;) = — f coth(A - uj) Z(X) dX . 
2m J 
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Then we obtain 



P„ ^ e/3^(«'+i«C) a = 0,±l. 



Observe that the function z has the following properties 

z{w + iQ — z{w — iC) = 1 — Z{w), 
z+{w) - z^{w) = Z{w), w£[-q,q\, 



(5.12) 



(5.13) 



where (11.14p is used in the first equation, and z± are limiting values of z from the upper and 
lower half planes in the second equation. 

Therefore, the leading non-oscillating part of the generating function can be written as 



G(0)(/3,m) =^(/3) •e^'"^r27r smh{2q) p{q) m\'^G'^[ 1,^^ e^^^°-^^^ 



(5.14) 



Here the constant Cq given by (|5.9|) . G(l, z) is the product of the Barnes functions (see (|3.27p ). 
and Ci reads 



C. = 1 / - dAdM + 2Z / ^ dA . 

2 J tanh(A — /i) J tanh(g — A) 



(5.15) 



The coefficient A{I3) is 



(e'3-l)2.det [/+A.f7W(^^^/) 



• det 



where 



e^^'("')[i^«(w - w') - K^{9i - w')] 



and 



?7gV(ii',^i' ) = - 



g-/3f(«)'-iC) _ g/3-/32(«)'+jC) 



(5.16) 



(5.17) 



(5.18) 



and the kernels Ug^ {w, w') and Ug^ {w, w') act on the contour T surrounding the interval [—q, q]. 
5.1.2 Leading oscillating behavior of (e^®™) 

Consider for exemple the function G^~^\p,m), the case G^^\(3,m) being completely analogous. 
Like in the previous case, the function (|4.72|) becomes linear at 7 = 1, and (f^^)' = l/27ri. 
Therefore the Fredholm determinant in (|4.78|) is equal to det[I + -^K] like for the non-oscillating 
term G^^\p,m) and the integral equation for p'{X) (15.41) remains the same. 
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As for the linear integral equation for z^^\X), it reads 

-q 

Using the expression (I4.46P of ^(A), one can write (I5.19P in the form 

q 

-q 

and the solution is 

z(+)(A) = Z(A)(^ + l) -1. (5.21) 

One should now notice that, everywhere but in the argument of F^^\ the function z^^\X) 
enters in (I4.78P only through the combination 

(^3 + 2^u)Z(A)_ 

which means that this part of the result can merely be obtained from the non-oscillating term 
via the replacement /3 ^ /? + 27ri. It is easy to see that it is enough to make the same shift 
P f3 + 27ri in the function : indeed, the Predholm determinants in this function contain 
the following products 



= "";;S" + ' + n " ^ > «=0'±i' (5-23) 

smh(t(; — q + laQ) smh(ty — Zs^p + laQ) 

s,p 

which, in the linearized limit {zs^p} = {As^p + es,p} become 

r^(+) \ sr^ / , • AX 1 sinhiw + q + iaO \ 

P^' ^ exp < > es^p coth(u; - A^^p + laQ + log — , . > 

smh{w — q + laC,) j 

^exp|y"coth(u;-A + iaC)(2^+^(A) + l)dA| = e.^P+^-^) K^+^O ^ (5.24) 

and one can similarly show that the same shift of /3 has to be done in the equation (15. Sp . 
The term G'^~\(3,m) can be considered in the same way, and we finally obtain 

GW(/3,m) = G(°)(/3±27ri,m). (5.25) 

It is worth mentioning at this point that the generating function (e^^™) is a polynomial of e^, 
which means that the exact result should be a 27ri-periodic function of (3. In the asymptotic 
formula this periodicity may of course be broken. We see, however, that the leading oscillating 
part of the asymptotics partly restores the original periodicity. It is therefore very possible that 
the more rapidly oscillating corrections to this formula can simply be obtained, at their leading 
order, by the shifts /3 ^ /? + 27rm, n G Z*, in G(°)(/?,m). 
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5.1.3 Final result and comments 



We eventually obtain 



GW(/5 + 27ria,m)[l + o(l)], (5.26) 



cr=0,± 



m which G(°)(/3, m), given by (|5.14p . corresponds to the leading non-oscillating term whereas 
27ri, m) are the leading oscillating ones. 
To conclude this section, we would like to stress once more that the result (|5.26p for the 
leading asymptotics of the generating function (e^®™) is formulated in terms of the solution of 
(11.14p . which is a /mear integral equation, whereas the integral equations describing the partial 
sums Gi,,,n{l) were non-linear. The linearization arises only after taking the last sum over n, 
in other words, only if we take into account the contributions of cycle integrals of all possible 
lengths. In the framework of our approach, the necessity to consider cycles of arbitrary lengths 
is quite natural. Indeed, the asymptotics of all cycle integrals have a rather common form 
independently on the length of the cycle. In particular it always contains a term which is linear 
over the distance m, and one might therefore expect that all cycle integrals eventually give a 
contribution of the same order to the asymptotics of (e^^™). It is worth mentioning however 
that, for cycle integrals of length £ > 1, this linear dependence on the distance shows itself only 
in the asymptotics, while it is explicit for the cycle integral of length 1=1: 



^^J^^H z sinh^(.-A) = / ^ H^o(A) + ^e] gi 



(5.27) 



Therefore, dealing with the series for (e^^™) in its initial form (|2.40p . we see only the explicit lin- 
ear dependence in m generated by each double pole at Zj = Xj (i.e. by the cycle integral of length 
i = 1), but we do not see the hidden linear dependence produced by all other cycle integrals. 
This may create the fallacious impression that the role of these double poles predominates. 

We specify here all these details since they explain the appearance of a non-linear integral 
dressing equation in the works [39l[38]. Exactly the same non-linear equation can be obtained 
from the series (|2.40l) after summing up the contributions of only the double poles. We would like 
to point out that this way leads to a misleading representation for the generating function (e^^""). 
In particular, using such representation, the authors of [39] concluded that the oscillations in 
the asymptotics of the two-point correlation function (crfa^^i) are exponentially suppressed at 
long distance, which is not true. 

5.2 Asymptotic expansion of (crfcr^„+i) 

Let us now consider the asymptotic behavior of the two-point correlation function {(^iC^m+i) 
which can be obtained via (11.41): 



2 



-4D + 1. (5.28) 

/3=0 
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Similarly as for (e^^™), we will consider separately the non-oscillating and oscillating parts. 
5.2.1 Leading non-oscillating behavior of (frfcr^^^) 

To obtain the leading non-oscillating term of the asymptotics of (cif cj^_,_]^), we should apply 
the differential operator 2D'^dp to the term G^^\(3,m) in (|5.1|1 . Evidently, in order to obtain 
non-vanishing contributions, one has to apply the derivatives over /? either to e^™^ or to the 
fractional power of m, setting simply /3 = in the remaining part. This gives us, for m — > oo. 



2Di 



df3'- 



m] 



/3=0 



2Z' 



A{/3 = 0). 



(5.29) 



Consider now the behavior of the factor at /3 — > 0. Let us take, for example, the deter- 
minant of the operator U^-^'^ . Due to the factor e^^^^'^ the kernel Ug^\w,w') has a cut on the 
interval [—q,q], therefore, for /? small enough, the action of the integral operator on the closed 
contour F can be reduced to 



(A)s 



-9.9] 



(5.30) 



where are the limiting values of U^'^^ from the upper (lower) half-planes. Using the equations 
(15.13P we obtain 



det 



det 



2TTi 



-<?.<?] 



where the operator in the r.h.s. acts on [—q, q\ and its kernel is 



(5.31) 



(5.32) 



Setting now /3 = we obtain 
1 



det 



det 



/3=0 



/ + —[K{w - w') - K{9i - w')] 
27r 



On the other hand. 



(5.33) 



det 



I + — [K{w - w') - K{9i - w')] 
2tt 



det 



I + —K{w -w') 
27r 



det [I-R{ei,w') 



det 



I + —K 

2-n 



R{ei,w)dw 



Z(9i)- det 



I + —K 

2-K 



(5.34) 



where R{w^ w') is the resolvent of the operator / -|- -^K. Thus 



det 



1 

2Tii' 



I + —Ui^\w,w') 



Z{ei) ■ det 



/3=0 



I+—K 
2-K 



(5.35) 
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One can prove similarly that 



det 



Z{e2)- det 



/3=0 



I + —K 
2tt 



Substituting these equations into (|5.16|) and using again (I5.13|) . we find that 
lim A(f3) = 1. 



(5.36) 



(5.37) 



Thus, taking into account (|5.28p and (I5.29p . we obtain the following leading non-oscillating 
asymptotic behavior of the correlation function: 

22:2 

■ ~ ~ ' (2Z? - 1)2 - + o 



l"m+l/ non-osc 



1 



m — > oo. 



(5.38) 



5.2.2 Leading oscillating behavior of {(yla^_^i) 

As we have seen, in order to obtain the leading oscillating term of the asymptotics of the 
correlation function, it is enough to shift /? by zb27ri in (|5.14|) . Observe that the constant A{l3) 
is proportional to (e^ — 1)^, therefore ^(±27ri) = ^'(zb27ri) = 0. Thus, in order to obtain a 
non-zero contribution, one should differentiate only the factor (e'^ — 1)^ when taking the second 
derivative with respect to j3, setting (3 = ±.2-Ki in the rest of (I5.14|) . 
Let us consider, for example, the case [5 ^ j3 + 27ri. It gives 



df3' 



m] 



A"{27ri) 



Ci—Co 



f3=0 



g2imp^ G2(l,Z)e 

[2tt sinh(2(7) p{q) m]^^^ 



(5.39) 



Here Pp is the Fermi momentum Pp = p{q) and 



2 det 



A"{2Tri) 



I + hK^i^^ ^Ol • det [/ + ^^0 



-. (5.40) 



|^g27ri5(6»i-HO _ g27ri2(ei-jC)j |^g-27rj5(6»2-jC) - g-27ri5 (02-^0] det^[/ + ^K] ' 

Recall that 0i^2 in (|5.40|) are arbitrary complex numbers. Let us set 6i = —q and 82 = 
Then it is easy to see that, at /3 = 27ri, 



2TTi 



(5.41) 



where f stays for Hermitian conjugation. Hence, the two determinants in the numerator of (|5.40p 
are complex conjugated. It is also easy to show using (|5.13p that 



^2mz{-q+iC) _ ^2mz{-q-iC,) 



^-2mz{q-iC) _ ^-2mz{q+iC) 



-4 sin^ t:Z ■ e2Ti[5(-g-iC)-5(g-iC)] _ 



^The ^-independence of the expression given in l|5.40p is proven in Appendix lA.3l This specific choice of 9j is 
the most convenient for the calculation of the Predholm determinants in l|5.40p in the vicinity of the free fermion 
point, see Appendix [B1 
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Combining (fOoD with the Barnes function in (09]l and defining A = -AG'^{l,Z) A!' {2m), we 
find 



A 



g7rj[2(g-iC)-2(-g-iC)] 



G(2,Z).det I + ^U^}l{w,w') 



-kZ ■ det [/ + ^K] 



(5.42) 



where G(2, Z) = G{2 + Z)G{2 - Z) and where we have used dOTD . 

One can transform similarly the expression for A!'{ — 2tti). The result in this case is still 
given by the same constant A (|5.42|) . Then, taking the second lattice derivative, we eventually 
obtain in the limit m ^ oo, 

[27r smh(2gj p{q) mj^^ 

The obtained amplitude of this leading oscillating term appears to be closely related to a 
special form factor of the operator . We will comment this relationship in the conclusion. 

5.2.3 Final result and comments 

Thus, at m ^ cxD, the two-point correlation function of the third components of spin in the 
external magnetic field behaves as 

(-f<+i)>.... = {^D - 1)^ - S + 8le^-^o sin^p, • . "^.^ ■ (5.44) 

^z^z y2,Ti smh(2gj p{q) mj^^ 



Note that, dependending on the value on A (and thus on Z), the second term of (I5.44p may be 
dominant compared to the third one or vice-versa. 

Let us discuss this result from the viewpoint of its comparison with known results and 
predictions. 

As already mentioned in the introduction, the asymptotic equation (I5.44p completely con- 
firms the predictions from Luttinger liquid and conformal field theory approaches. In the limit 
of free fermions C = f (see Appendix [B|), this result also agrees with the known answer 

{afa^m+i)A=o = (2^ - 1)' - - cos(2mp^)]. (5.45) 

In the general case, there exists no prediction concerning the value of the amplitude of the 
oscillating term, except in the zero magnetic field limit (see [Ml ESI [66]). However, in the 
framework of our approach, the non-zero magnetic field plays the role of a certain regularization 
which ensures the finiteness of all the constants (Co, Ci, A, p{q) sinh2g) entering this amplitude: 
at /i = all these constants become divergent. It is nevertheless possible to show that the total 
combination of these constants remains finite. As we have mentioned already the corresponding 
proof is highly non-trivial, therefore we do not present it here. We would just like to mention 
that the most complicated part of the calculations is related to the extraction of the divergent 
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part from the Predholm determinant det[I + -^U^g {w , w')]: one can show that it coincides with 
the divergent part of det[/ — ^K]. 

Due to the difficulties that arise when taking this hmit /i — > in our result, the complete 
comparison of (|5.44|) with the results of [Ml EH [66] still remains an unsolved problem. Note 
nevertheless that, in the limit of free fermions, the formula (I5.44P holds for arbitrary magnetic 
field including the case h = 0. We succeeded moreover to compute explicitly the limiting value 
at /i = of the amplitude of the oscillating term in the vicinity of free fermions up to the second 
order in e = - — i. Our result, 



2(-l)'" 
{(^tcrm+i)osc = n exp 



4e(log2 + C) + 8e2(2f -l-C-log2) +0(e^), (5.46) 



6 

C being the Euler constant, coincides at this order in e with the one of [Ml [65l [66] . 

6 Quantum non-linear Schrodinger equation 

In this section we briefly explain how to apply the method described above to another model, 
the quantum one-dimensional Bose gas (or quantum nonlinear Schrodinger equation model). 

Recall that the starting point of our derivation is the master equation (|2.ip . We have 
stressed already that such integral representation exists not only for the XXZ chain, but also 
for other models possessing the six-vertex i?-matrix [IS] and, in particular, for the system of 
one-dimensional interacting bosons on some finite interval [0,L]. 

The Hamiltonian of this model is given by 

H=[ (d^.^'^d^^ + c^''^^^^ -h^^^)dx. (6.1) 



Jo ^ ' 

Here ^ and ^'^ are Bose-fields possessing canonical equal-time commutation relations, c is a 
coupling constant and h a chemical potential. For c > and /i > 0, the ground state of the 
model goes to a Dirac sea in the thermodynamic limit and can be described by a set of integral 
equations similar to those of the XXZ chain. 
The analog of the Qm operator is 

Q,= r <^\z)^{z)dz, (6.2) 







and the function (e^^"=) is a generating function for the correlation function of the densities 
j(x) = ^'^{x)^{x) via 

1 ^2 92 



(6.3) 

/3=0 



The method to compute the asymptotics of (16. 3|) in the case of the one-dimensional Bose 
gas literally repeats the steps we have described in details in the case of the XXZ chain. Let us 
just point a few peculiarities. 
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The master equation for (e^®^) has the form (|2.1|) . but the functions a(A), d{X), and /(A) 
(see (1231, (I23|)) are now 

a(A) = e-'^, (i(A)=e'^, /(A) = e"^^^. (6.4) 



In the remaining part of the master equation, all hyperbolic functions should be replaced by 
rational ones, i.e. sinh(A — z ± i() ^ \ — z ± iC, etc. After this operation one should set ( = c. 
The transformations of the master equation described in Sections 12.21 and 12.31 are the same for 
the one-dimensional bosons. 

One of the difference between the quantum one-dimensional Bose gas and the XXZ model 
is due to the choice of notation^]. In particular, the integral equation for the dressed charge in 
the Bose gas model coincides with its XXZ analog up to the sign of the integral operator: 
g 

Z{X) - ^ y K(A - /z) Z(/.) d// = 1 , K{X) = -j^^ . (6.5) 



One should also pay attention to the fact that the function /(A) (|6.4I) can be analytically con- 
tinued to the lower half-plane, while its analog (|2.5p can be continued to the upper half plane 
(in the case of ivr-periodic hyperbolic functions, we can of course always consider the domain 
[Q'AI < f ). This difference leads to some differences of signs in several formulas, but the common 
strategy of the derivation remains just the same as in the case of the XXZ chain. 
The final answer for the leading asymptotic behavior of ()6.3I1 is 

(i(.)i(0)>„,„ -r,-'-^^ 2A,l . '^"^^ . 

Here, just as in the previous case, Z = Z(ibg), Z{X) satisfies the integral equation (I6.5|l and q, 
the value of the spectral parameter at the Fermi boundary, depends on c and h (see e.g. [59]). 
In the model of the one-dimensional Bose gas the spectral density p{X) is related to the dressed 
charge by p = ^ . The Fermi momentum Pp and the average density D are given by p.lSjl . 
The expressions for the constants Co and Ci are quite similar to (|5.9I) . ()5.15|1 . 

' Z(A)Z(/z) 
{X — fi — icY 



Co= / TV^^^z^dAd//, (6.7) 



= 1 / ''^'^ '^^\ - '^'^ ''^^^ dA d, + 2Z / ^ - ^1^) dA . (6. 
2 J X- II J q-X 

-q -q 

The formula for the constant A coincides also basically with (I5.42p : 



A 



_,m\z{q-tc)-z(—q-ic)\ ^ 



2 



(6.9) 



^One could avoid this difference by replacing ^ tt — ^ in the original formulas. 
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Here 

9 



z{w) = — I -^^dX, (6.10) 

^ ^ 27ri X-w ' ^ ' 



and 



^-"^^^ = %2...(.+l) _ e2..i(.-.c) [^(^ - - ^('? + ^0] . (6.11) 

It is not difficult to check that the result (I6.6|l has the correct free fermion limit (c = oo). 
One can also show that it reproduces at first order in the asymptotic behavior obtained in 

EHl. 



Conclusion 

We have described a method to compute the asymptotic behavior of correlation functions of 
quantum integrable systems. Our study is based on the multiple integral representation (12. 
called the master equation. Such master equation can be obtained via the algebraic Bethe 
ansatz for a rather wide class of integrable models and different correlation functions. In the 
present paper, we have considered only one specific correlation function for the XXZ chain and 
one-dimensional bosons. To conclude, we would like to discuss possible further developments. 

First of all, it should be possible to apply this method to compute the asymptotic behavior of 
other correlation functions, such as {(Tfa^_^_i) for the XXZ chain or (^'^^(x)^'(O)) for the Bose gas. 
It would also be interesting to apply this approach to compute the long-distance correlations in 
the massive phase A > 1 of the XXZ chain. In that case there is no power law corrections in the 
asymptotic behavior of the cycles. This should lead to an exponential decay of the correlations. 

The case of the XXX model is also of special interest. For non-zero magnetic field the result 
(|5.44p apparently remains valid. However, at /i = 0, one can expect logarithmic contributions 
to the asymptotic expression OlS]. It would be interesting to see how such terms can appear 
in the framework of our approach. However, the case of zero magnetic field appears to be the 
most complicated one, even for the XXZ chain. Therefore before studying the XXX model it 
would be desirable to understand completely what happens for the XXZ chain at h = 0, and in 
particular to compute explicitely the (finite) limiting value for the amplitude of the oscillating 
term. 

We hope that our approach can be used for the asymptotic analysis of the temperature 
dependent correlation functions. At least in the case of the one-dimensional Bose gas, for which 
there is no bound states in the spectrum of the Hamiltonian for c > 0, the generalization of 
our method looks quite straightforward. For T > 0, the integrals over the spectral parameters 
should be taken over the whole real axis. This leads to the absence of power law corrections in 
the asymptotic behavior of the cycle integrals, just as in the massive regime of the XXZ chain. 
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As a result, the asymptotic behavior of the correlation functions should decrease exponentially 
with the distance. 

We would like finally to discuss the nice relation that exists between the amplitude of the 
leading oscillating term of the asymptotics (|5.44p and a special form factor of the operator , 
namely the matrix element 

\2^) ■ ||VXM)||-!l^({A})i|' ^'-^'^ 

where is the A^-particle ground state and is an excited state containing one 

particle and one hole at the different boundaries of the Fermi sphere, for instance, Ap = q, 
\h = —q- We recall that in (|6.12|) M denotes the length of the chain. Using the determinant 
representations of form factors obtained in [52], one can show [U] that, in the thermodynamic 
limit |-F(jP is equal to, 

lim = AA sin^p^ (f^-Co ]^2tt sinh(2g) p{q)]''^^\ (6.13) 

Af^oo 

Thus, the amplitude of the leading oscillating term is equal to the square of the norm of the 
form factor ()6.12p : 



K^m+l),eading = (2^ - 1) -^^+2|i^<x| " ^2^2 • (6.14) 

The coefficient 2 is due to the fact that there are two such form factors. This observation 
perfectly agrees with the conformal field theory approach. Namely, if we consider the two-point 
correlation function as a sum of form factors, then the main contribution to the asymptotics 
comes from the terms corresponding to the excitations in a vicinity of the Fermi boundaries. 
Hereby, the leading oscillating term in the asymptotics is produced by the excited states having 
particles and holes on the different sides of the Fermi zone. Thus there is an opportunity to 
improve drastically the form factor approach. This, in turn, would open a path towards the 
computation of dynamical correlation functions. 
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A Determinant representations 
A.l Proof of Proposition [2TT] 

If the parameters Ai,...,AAr satisfy the system of Bethe equations (12.21) . then d{Xj) can be 
expressed in terms of a{Xj), and the determinant of can be written as 

N N 

detNn,{{z},{X}\{z}) =lla{Xj)- J] sinh(z, - - <) ' det^ M«({z}|{A}), (A.l) 



a,b=l 



where 



(M«) = t{zk, Xj) + Kt{Xj,Zk) ■ V+{Xj) Vl\Xj) 
In order to prove (|2.10l) . we shall prove that 

N 



detjy = det 



N 



1 



sinh(Aj - Zk) 



n 



1 - K 



Consider the following transformation of the l.h.s. of (IA.3P 

detiv(M«A) 



• det 



N 



det jv 



det AT A 



with 



N 



n sinh(2j - Xa) ( coth(zj - A^) for k^N, 



n sinh(zj - Za) I ^ 



for k = N. 



a=l 



The determinant of A can be easily computed. Indeed, 



N 

n sinh(za - Xb) 
det N A = ^— • det m 

N 

n sinh(2:a - z^) 
a,fe=l 



/ coth(zi - Xk) 
\ COth(zAr - Xk) 



(A.2) 



^^k + u\l\e)\. (A.3) 



(A.4) 



(A.5) 



(A.6) 



Subtracting the last line from all others, we reduce the remaining determinant to a Cauchy 
determinant, which gives 



N 



det TV ^ 



n sinh(za - Aat) at 
a=i -rr smh(Aa - X^) 

n sinh(A, - Xn) ^ ' 

a=l 



(A.7) 
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The product of matrices Mi^A can also be explicitly calculated 



•jk 



Si 



N 

n sinh(Aj - Xa) 

a=l 

a-+3 



Y\ sinh(Aj 



K^{Xj-Xk)V+{X,), k<N, (Ai 



a=l 



and 



iM^A),N = il-K)V+iXj), 



N. 



(A.9) 



Here the function Kk{X) is given by (|2.15|) . Let us explain how the formulas (jA.Sp . (|A.9p were 
obtained. For instance, in order to obtain (lA.Sj) . one has to calculate the sums 



N 

E 



—i sm( coth(2;^ — Xk) 
^ sinh(z£ — Xj) sinh(z^ — Xj ± i() 



N 

n sinh(2;£ - Xa) 

a=l 

N 

n sinh(2;£ - Za) 

a=l 



(A.IO) 



Consider an auxiliary contour integral 



-i sinC coth(?i; — A^) dw 



liri J sinh(tt; — Xj) smh{w — Xj ± iQ sinh(t(; — Za) ' 



N 

n 



sinh(?i; — Xa) 



(A.ll) 



where the integration is taken over the boundary of a horizontal strip of width in. On the 
one hand, the integral vanishes I-t = as the integrand is an ivr-periodic function that is 
exponentially decreasing for ^{w) — > ±00. On the other hand, it is equal to the sum of the 
residues inside of the integration contour. The sum of the residues at the poles w = Za gives 
exactly G^^. Taking into account the contribution of the additional poles at w = Xj =F and 
w = Xj (the last one only exists for j = k), we arrive at the following identity 

N 



n sinh(Aa - Xj) 

a=l 

G% ± Sjk ^ ± coth(A,- -XkT iC) ^t(Aj) = 0. 



(A.12) 



n sinh(za - Xj) 



a=l 



We have thus computed and in this way proved the formula (lA.Sjl . The result (|A.9p can be 
obtained by a similar method. 

The following transformations are trivial. We can extract the factor V+{Xj) from each line 
of (M«;A)jfc, making the elements of the last column equal to 1 — k. After this, subtracting 
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the last line from all the others and extracting the coefficients in front of 6jk we obtain a new 
representation for det 



1 - K 



N 

X det AT 



V_(A„) 
1 



sinh(Aj — Zk) 



■ det 



N 



5jk + u\^^{\N)\. (A. 13) 



Thus, we have reproduced the r.h.s. of (|A.3p up to replacement 6 — > \n 

Now consider det 
explicit form). After a simple similarity transformation we obtain 



5jk + U^k* {9) appearing in the r.h.s. of equation (|A.3p (see (|2.13p for its 



N 

a=i Ki^{Xj - Xk) - K^{6 — Xk) 



n smh(Aa - Afc) + ^ ^ 

a=l 



(A.14) 



Let us multiply the first — 1 columns by the coefficients Sk 

V^\Xn) - kVi\Xn) 
and add them to the A^-th column. Then the last column of the determinant becomes 



(A.15) 



N-l 



k=l 



v:;:\e) - kVi\9) 

V^\Xn) - ^^VI^XnY 



(A.16) 



(the method of calculation is quite similar to one described in the formulas (IA.10j) - (|A.12h ). Now 
it is enough to subtract the last line of the obtained matrix from all others, and we arrive at the 
following identity 



det 



N 



V^\Xn) - kV-\Xn) 
In other words, it means that the combination 



det AT 



(A.17) 



det 



N 



(A.18) 



does not depend on 9. Thus, substituting (lA.lTj) into (|A.13p . we obtain the r.h.s. of (IA.3j) . 
The representation (|2.1ip for can be proved by using the identity 

detTvM^ = det7vM«, where {M^) = t{zk, Xj) + ^^t{Xj, Zk) ■ V+izk)Vl\zk), (A.19) 
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established in [35] (see Appendix B of this paper). It follows from (|A.19|) that 



N 



detjvM«({4,{A}) =^'^ n 



sinh(Aa - Zb- i() 



•detivM^-i({A},{z}). 



(A.20) 



^-^^^ sinh(zf, - Xa- iC) 
Then the representation (|2.11|) follows from (I2.10p after the replacements z -s-^ A and k 



A. 2 Fredholm determinant representation 



Consider the Fredholm determinants det 



I + ^^U^g^''\w,w') (see (12^6)1 ). We have 



log det 



E 



' k J (27ri) 
^ r 



U',^''>{wuW2) ■ ■ ■ U^,^''\wk,wi). (A.21) 



Computing the multiple integrals by the residues at Wj = (respectively at Wk 
i = 1, . . . , N , we obtain 



log det 



k=l 



(-1) 



fc+1 



N 



f A z) 

what is exactly the expansion for log det at Sjk + Ujj^' {6) 



where 



(A.22) 



A. 3 Some identities for Fredholm determinants 

Here we prove ^-independence of Fredholm determinants considered in Sections [2] and El 
Consider an integral operator I + V acting on some contour V . 

Proposition A.l. Let 



g{w) = h{w) + I V{w,w') h{w') dw' 



Then, for arbitrary wq € C, the following identity holds 



det \I + V{w,w')] =^^.det 
h{wo) 



I + V{w,w')-4^V{wo,w') 



(A.23) 



(A.24) 



Proof — Suppose that the resolvent R{w,w') of the operator I + V exists (if not, then we 
can consider some regularization Ve{w,w') of the original kernel, such that the corresponding 
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resolvent exists). Then 



h{wo) 
1 



det 



I + V{w, w') - ^P^V{wo, w') 



det [I + V{w,w')] ■ det 



9{w) 
9 (wo) 



R{wo,w'] 



■det [I + V{w,w')] (giwo)- / R{wo,w')g{w') dw' 



= det[l + V{w,w')] , (A.25) 

which ends the proof. □ 

If the kernel V{w,w') is analytic within some domain V containing F, then the identity 
(1X241) holds for woeV. 

Let now T be the contour shown on Fig. 1 and h{'w) an ivr-periodic function that is holomor- 
phic outside of T and bounded at tt; — > ±00. Then 

/dw' 
\KJw-w')-KJe-w')]h(w') = Kh(e-iC)-h(e + iC)-Kh(w-iC) + h(w + iC), (A.26) 
2m 

r 

where the integral has been computed by the residues lying outside of the contour F. Compose 
a kernel u'^\w,w') as 

h{w) 



■ [K^iw-w')-K,{e-w')]. 



2-111 [Kh{w — iC,) — h{w + iC)] 
Let us apply (|A.24ll to the operator I + U^g'^\w,w') with 6* € F. We have 



g{w) = h{w) + 6 Uf\w,w') h{w') dw' = h{uj) 



Kh{e - iC) - h{e + iCi 

Kh{w — iC) — h{w + iQ 
r 

Substituting this into (IA.24P we immediately arrive at 



(A.27) 



(A.28) 



det 


/4 


-^^U^'\W,W') 


Kh{e - 


iC)-h{e + iC) 



det 



I+^,uL':^{nj,w') 



(A.29) 



for arbitrary wq such that |9(tyo — < C, w G T . 
In particular one can take 



N 



smh(w - \a) 

a=l 



(A.30) 



In the first case we obtain the identity for the operator U^^\ in the second for U^'^\ Similarly 
one can prove ^-independence of the kernels JJ^^^ and U'^^^ . 
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B Free fermions 

Consider the determinant of the matrix Sjk + Uj^\6) in the limit of free fermions C = f- Then 
V-^-{f^) = and the entries {9) are equal to 



N 



J] tanh(za - Aj) 



^ [tanh(Aj - Afc) - tanh(6' - Afc)] , 

n tanh(Aa - \j) 

a=l 



(B.n 



The straightforward calculation of the corresponding determinant causes serious difficulties. 
However, there exists a way to avoid these problems. Due to (IA.3I) the determinant of + 
Uj^\6) can be expressed in terms of detM^. One has, for the limit of free fermions. 



det 



N 



det 



N 



1-K 



(B.2) 



_sinh{2fe-Aj)_ 

On the other hand, the matrix (|A.2p becomes the Cauchy matrix at C = f , 
2(1 - k) 



vr 



for C = — • 

^ 2 



■^'^ sinh2(2;fc — Aj) ' 

Thus, both determinants in the r.h.s. of (IB.2P are explicitly computable and we arrive at 

N 



(B.3) 



det 



N 



N 



Similarly one has 



det 



N 



N 



cosh(0 — 


Za) 


n cosh(2;a - Zh) cosh(A„ 

a>b 


-Xb) 


cosh(0 — 


Xa) 


N 

n COsh(Za - Afe) 
a,b=l 




cosh(0 — 


Aa) 


N 

n COsh(Za - Zb) cosh(Aa 
a>b 


-Xb) 


cosh(6' — 


Za) 


N 

n cosh(2;a - Xb) 





vr 



C = 2 • (B-4) 



(B.5) 



a,b=l 



It is easy to see from the obtained result and the definition (|2.2Qp of Wn that, in the free fermion 
limit. 



Wn 



{A} 

{z} 



1. 



(B.6) 
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Let us now compute the constant A in the hmit of free fermions. For ^ = | we have Z{X) = 1, 
hence, 

U^g{w, w') = — tanh(t(; — q) coth(w + q) [tanh(t(; — w') + tanh(g + w')\ . (B.7) 

The kernel (IB.7P is holomorphic in a vicinity of the interval [—q-,q\. We draw the reader's 
attention on the fact that this property only holds due to the special choice of 6i = —q. Otherwise 
the kernel Ug^\w,w') would have a simple pole at w = —q. 

Since the integral operator {w,w') acts on the closed contour surrounding the interval 
[—q, q] we obtain that 



det 
hence, 



1, (B.8) 



^ = J_g2Ti[2(g-j|)-2(-g-if )]_ ^g_g^ 

It is also easy to see that for free fermions Ci = (see (I5.15P ). and 

2Tri[z{q-f)-~z{-q-f)]-Co = 0. (B.IO) 
Thus, taking into account that p{X) = ^cosh(2A) for ^ = f , we obtain 



, z z \ _ f sinpp Y 2cos(2mpj,) 
It remains to observe that 

9 



-Q 



(B.ll) 



/d A 
^^^^^ = arctan (e^^) - arctan (e"^^) . (B.12) 



Prom this we find that sinp^ = tanh(2(7), and hence 
,zz N 2cos(2mp^) 

C The Lagrange series and its generalizations 

In this appendix we consider several generalizations of the Lagrange series (see e.g. [81j) used in 
Section 14.31 We give the detailed proof in the standard (scalar) case. The generalizations then 
become quite evident. 
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C.l Scalar case 

Let us consider the series 



°° 1 

n=0 



(C.l) 

e=0 



where F(e) and 4>{e) are holomorphic for |e[ < tq. 

Proposition C.l. // there exists r < such that \4>{e)\ < r for \e\ = r, then the series I C.l]) 
is absolutely convergent and its sum is given by 

_ Fjz) 

where z is the root of the equation 

z-(j){z) = (C.3) 
such that \z\ < r. 

Proof — Replacing the n}^ derivative by a Cauchy integral, we obtain 

~ ^ ^ (C.4, 



2vri , 

n=0 I "f 

a; =r 



Since |</>(w)j < r = the obtained series is absolutely convergent, and we arrive at 

Go = ^ / -^^,dz. (C.5) 

ZTTl J CO — (f)[uj) 
a;|=r 

Due to Rouche's Theorem, the equation u; — 4>{<jj) = has exactly one simple zero uj = z within 
the circle |a;| < r. Taking the residue in this point we obtain the statement of the proposition. □ 



C.2 Matrix case 

The result obtained for the series (IC.lj) can easily be generalized to the case of several variables. 
Consider a multiple series of the form 

where F({e}) and 4>j{{e}) are holomorphic for \ek\ < r^^\ k = I, . . . , N . 



68 



Proposition C.2. If there exist rj < Vj^^ such that < rj for \zj\ = rj, then the series 

lie. 6\) is absolutely convergent and its sum is given by 



G 



N 



(C.7) 



detAT Sjk 

where Zj are the roots of the system 

zj - (t>i ({4) = 0, (C.8) 
and detAT Sjk is the Jacobian of the system II C. 

Sjk = S,k - i-Jj {{z}) . (C.9) 

The proof is completely analogous to the one of Proposition IC.ll It uses the analog of 
Rouche's Theorem for several complex variables (see e.g. |i3j|). 

We consider now special cases of the series (IC.6|1 . Namely, let /zi, . . . , //at be a set of complex 
parameters. Let 



where 9{X,fi) is some smooth function. Let also 

/ N N \ 

\ a=l a,b=l / 



where g^^\g^'^^ ■ ■ ■ are smooth functions. Then the obtained result takes the form 



Gn 



N N 
F( g'^^KfJ-a) Za, E a'^^KfJ-a, fJ-b) Za Zb] 

.a=l a, 6=1 



det 



N 



^jk - Oit^k,lJ'j) f'(^ Y.^Za9{fla,l^j)^ 



where 



/ (^Ylza9{fia,f^j)^ = 0. 



(C.IO) 



(C.ll) 



(C.12) 



(C.13) 



C.3 Continuous case 

Let us now consider a series of multiple integrals 

q 



^ = E /d"A n ^ n / E ^(A- A.) • ^ E 9''' (A.) e.; 

n=0 i=l ^ j=l \a=l J \a=l 



(C.14) 



e,=0 
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where for brevity we have written exphcitly only the first argument of the function F. Consider 
a discreet analog of these multiple integrals, 



n=0 ■ Ai,...,A„eA i=l ^3=1 \a=\ / \a=l 



Here each \j independently runs through a finite set of values A = {/ii, . . . ,/iAr}, where /Ufc+i — 
\ik = A. Obviously G'(A) ^ G at A ^ 0. 

We now reorganize the inner sums over the lattice into sums over the number of A's equal 
to some point of the lattice. Namely, we sum up with respect to all the possible numbers Sj 
such that there are Sj A's equal to the lattice point /ij. Obviously there exist ways of 

realizing such a configuration. It is also evident that, for a configuration such that, for each j, 
there are Sj A's equal to /ij, then for any function one has 



d 



\a=\ 



N 



N 



, (C.15) 



Hence, 



N 



o(A)=EE n 



n=0 Si>0 j=l ] 



N 



N 



N 



\a=l 



\a=l 



T,Si=n 

Changing the order of summation and re-scaling ej — > Aej, we obtain the series (|C.6P with 
and F in the form (jClOj) . (IC.lip . in which the functions and 5^*) should be replaced by 
A9{fia,fJ.j) and Ag^'^^fia)- Then the continuous limit A ^ is trivial, and we finally obtain 



G 



F{ /5«(^)z(^)d^; 

-Q 



det 



(C.17) 



where the function z{fi) satisfies the integral equation 



z{t^) = f\ y 0(A,^)z(A)dA 1 . 



(C.18) 
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C.4 Multiple series of multiple integrals 

Consider now a multiple series of the form 



E nA/nn 



dA 



9e 



n £s 



n n E E ^( V' ^^.p) • ^ ( E E 5^'^ (^*.'^) • • • 

s=lp=l \t=l a=l / 



,t=l a=l 



. (C.19) 



Just like in the previous example we can take the lattice approximation of the integrals. Then 

it is easy to see that 



Gi...„(A) 



(C.20) 



where 



/ n N \ 
:,a - /s I A ^ ^ Zs^b Oif^b, ^^a) \ = 0, 
\ s=l h=l / 



(C.21) 



and deinN Sjk is the Jacobian of the system (|C.2ip . Taking the sum over s in (|C.21j) . we obtain 



TV 



Za - /s„ A ^ e{llb, ^a) = 0, where Za = ^ Zs,a, /e„ = E •^'^ 



b=l 



s=l 



.8 = 1 



(C.22) 



It is also clear that detmN Sjk = det]\f Sjk, where 

Sjk = Sjk - A 9{iik, IJ-j) /e„ ^A ^ Zb 9{nb, /Xj) j 
is the Jacobian of the system (|C.22jl . Thus, in the continuous limit, we have 



(C.23) 



Gl,,,n 



F[ /5«(^)zW(/i)d/i; ... 

\-q 



det 



5( A - ^) - A) J / e{u, A) {u) du 



(C.24) 



where /s„ = YTs=i fk, and 



^^"^(/^)=/eJ y ^(A,^)z(")(A)dA 



(C.25) 
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